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Abstract 

Background: Low Emission Zones (LEZs) are areas where the most polluting vehicles are restricted from entering. The 
effectiveness of LEZs to lower ambient exposures is under debate. This study focused on LEZs that restricted cars of Euro 1 
standard without appropriate retrofitting systems from entering and estimated LEZ effects on NO2, NO, and NOx { = N02+ 
NO). 

Mef/jocfe; Continuous half-hour and diffuse sampler 4-week average NO2, NO, and NOx concentrations measured inside and 
outside LEZs in 17 German cities of 6 federal states (2005-2009) were analysed as matched quadruplets (two pairs of 
simultaneously measured index values inside LEZ and reference values outside LEZ, one pair measured before and one after 
introducing LEZs with time differences that equal multiples of 364 days) by multiple linear and log-linear fixed-effects 
regression modelling (covariables: e.g., wind velocity, amount of precipitation, height of inversion base, school holidays, 
truck-free periods). Additionally, the continuous half-hour data was collapsed into 4-week averages and pooled with the 
diffuse sampler data to perform joint analysis. 

Resu/ts: More than 3,000,000 quadruplets of continuous measurements (half-hour averages) were identified at 38 index and 
45 reference stations. Pooling with diffuse sampler data from 15 index and 10 reference stations lead to more than 4,000 
quadruplets for joint analyses of 4-week averages. Mean LEZ effects on NO2, NO, and NOx concentrations (reductions) were 
estimated to be at most —2 \ig/m^ (or —4%). The 4-week averages of NO2 concentrations at index stations after LEZ 
introduction were 55 |xg/m^ (median and mean values) or 82 \ig/m^ {95th percentile). 

Conclusions: This is the first study investigating comprehensively the effectiveness of LEZs to reduce NO2, NO, and NOx 
concentrations controlling for most relevant potential confounders. Our analyses indicate that there is a statistically 
significant, but rather small reduction of NO2, NO, and NOx concentrations associated with LEZs. 
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introduction 

Low Emission Zones (LEZs) are areas or roads where the most 
polluting vehicles are restricted from entering. They are currently 
introduced in 13 European countries [1]. In Europe, vehicle 
emissions are classified by the so-called "Euro Standards" with a 
current range from Euro 1 to Euro 6 regarding the technical 



features of the vehicles which are fixed in several EU-Directives for 
passenger cars and heavy-duty trucks [e.g., 2]. Basically, this 
means that vehicles are restricted in relation to their Euro emission 
level. The configuration of LEZs is extremely different and 
heterogeneous in Europe, for example in Italy, where the entry 
standards, the subsistent regulations and the daily duration of LEZ 
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conditions differ substantially from town to town. However, most 
LEZs in Europe operate 24 hours a day, 365 days a year [see 1]. 

One of the most developed applications is found in Germany. 
Low emission zones have been introduced in Germany since 2008 
in different stages, resulting in meanwhile 48 LEZs with 
restrictions for pollutant groups 2 or 3 in 11 Federal states by 
the end of January 2014 [3]. In this study we analysed the effect of 
introducing the "LEZ of poUutant group 1" which restricts from 
entering Diesel cars of an European emission standard below Euro 
2 without particulate reduction system and gasoline cars of an 
European emission standard below Euro 1 without appropriate 
exhaust gas catalytic converters [3]. 

Traffic emissions are considered to be a relevant source of air 
pollution [4] and LEZs are believed to be the most effective 
measure that cities can take to reduce vehicle-induced air pollution 
problems in their area [5-7]. The emissions that are aimed to be 
reduced by LEZs are mainly fine particles like PMIO or smaller 
[8—12]. The effe( tivc'n(;ss of LEZs to reduce traffic-related 
exposures is stiU under debate [13] and there is an open discussion 
in the public about the "outcome" and cost-benefit ratio of LEZs 
[14—16]. Most of the published information refers to particulate 
matter. 

Additionally, nitrogen dioxide is discussed to be a major traffic- 
related poUutant as well as an epidemiologic marker of air quality 
and related adverse health effects [17-21]. On the other hand, a 
systematic literature review showed only moderate evidence for 
adverse health effects at a long-term exposure below an annual 
mean of 40 |ig7m'^ N02 [22]. 

According to EU rules [23,24] limits were additionally imposed 
for N02 and are enforced in Germany since 2010: 200 Hg/m^ as 
an 1 hour average (acceptable: 18 excursions/year) and 40 Hg/m^ 
as an annual average. Values were and are in excess: about 69% of 
all stations near to traffic showed annual averages higher than 
40 |J,g/m'^ in Germany [7,25]. This non-compliance is not 
restricted to Germany but the European limit value for N02 is 
exceeded in many European cities [26-28]. The LEZ concept was 
extended and it was assumed that LEZs are an effective measure 
not only to lower PMIO dust levels but also to reduce N02 
C(m[:entrations [6,29]. There are indications that LEZs may 
indeed reduce NOx concentrations effectively [30-32], but ozone 
has to be considered a confounder in N02 measurements [e. g. 
33], and the gases NO and N02 rapidly interconvert, too [34]. 
Furthermore, national emission ceilings were defined for NOx, 
i.e., the sum of N02 and NO [35]. Thus, there is interest in the 
impact of LEZs on concentrations of NO and NOx also [36] . 

However, a scientific proof of the LEZ concept targeting at 
NOj, NO, and NO^ is stiU missing. In order to test the views of 
legislators and researchers that LEZs are effective measures to 
reduce nitrogen oxide concentrations [29,37], this study focused 
on the potential effects of LEZs on ambient concentrations of 
NO2, NO, and NO, in LEZ areas of 1 7 German cities. 

We reported on the effect of LEZs on PMio concentrations 
elsewhere [32]. 

Methods 

Target parameters 

The aim of the study was to analyse the effectiveness of German 
LEZs (as many as eligible) to lower NOg, NO, and NO,, ( = NO2+ 
NO) concentrations. The first analysis series of NOg, NO, and 
NOx were based on continuous half-hour measurement data of 
NO2 and NO. Second, measurement data for NO2 and NO 
concentrations collected by diffuse samplers and determined over 
longer sampling periods were available. These data were allocated 



to 4-week periods. Third, we collapsed the half-hour measurement 
data to four-week averages and pooled these collapsed continuous 
data and the diffuse sampler data to perform joint analyses over 4- 
week periods. The original NO2 and NO measurements were 
performed by the Environmental State Institutions in Germany 
(Landesumweltamter). A federal data base [38] reports on the 
applied measurement procedures. 

Measuring procedure 

Two measuring procedures were applied: continuous measure- 
ment devices (chemiluminescence), data stored as half-hour 
averages and diffuse samplers (Palmes tubes, chromatography), 
data stored as long-term averages over weeks. The chemilumi- 
nescence method relies on the reaction of NO with 03: NO+ 
03^N02*+02. Chemiluminescence is generated in the range of 
600 nm to 3,000 nm when the excited molecules return to the 
ground state. The light intensity is proportional to the concentra- 
tion of NO molecules. A deoxidation comx-rter is used to reduce 
N02 to NO. Thus, the N02 concentration is determined as the 
difference between the NOx concentration measured when the 
sample gas is directed through a deoxidation converter and the 
NO concentration measured when the gas is not run through the 
com erter. The diffuse samplers were Palmes type tubes modified 
with a glass frit as turbulence barrier. In these passive samplers 
molecules diffuse because of a concentration gradient through an 
intake opening with a defined cross-section along a fixed diffusion 
path to a sampling medium by which they are adsorbed. This 
process is described by Pick's first diffusion lnw. The chemical 
analysis is done by chromatography. [More details on both 
methods may be found in 39,40-43] . 

Period of investigation 

The period of investigation was from 2005 until the end of 2009 
(3 1 December 2009), starting at least from the introduction of the 
individual LEZ minus the length of the respective LEZ phase (or 
earlier if restrictions of truck traffic were enforced before the 
introduction of the LEZ). 

Low Emission Zones 

There were 34 German active LEZs until the end of 2009 and 
774 monitoring stations in use. With introduction of these LEZs, 
as a main effect, only those diesel vehicles with an exhaust 

emission standard better than Euro 1 (with sticker) were allowed to 
enter the zone. In principle, the German "LEZ of poUutant group 
1" restricts from entering 

— Diesel passenger cars, trucks and buses of an European 
emission standard below Euro 2 without particulate reduction 
system, and 

— Gasoline passenger cars, trucks and buses of an European 
emission standard below Euro 1 without appropriate exhaust 
gas catalytic converters. 

Local authorities can set up exception permits especiaUy for 
fight duty vehicles, trucks and buses due to local necessities [3] . 

According to protocol LEZs were included into the study if and 
only if 

— monitoring stations existed, that operated before and after the 
LEZ introduction and measured inside the LEZ area {index 
stations) and 

— monitoring stations existed, that operated before and after the 
LEZ introduction and measured outside the LEZ area - in a 
circle around the centre with a radius of about 25 km - and if 
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outside the city area, than in no other LEZ {reference stations) 
and 

- these monitoring stations measured NO2 or NO (continuous 
measurements or diffuse samplers). 

(For the terminology and the use of index and reference values 
in comparisons if exposures levels see Rothman et al. [44]) 

Seventeen cities with LEZs in 6 German Federal states could be 
included into the study (Baden-Wiirttemberg: Herrenberg, 
Ludwigsburg, Mannheim, Reutlingen, Stuttgart, Tubingen; Ba- 
varia: Augsburg, Munich; Berlin: Berlin; Hesse: Frankfurt; Lower 
Saxony: Hannover; North Rhine-Westphalia: Dortmund, Duis- 
burg, Dtlsseldorf, Essen, Cologne, Wuppertal). Figure 1 shows all 
active 34 German active LEZs in December 2009 and the 17 
LEZs included for study. File S 1 entails maps of all LEZs eligible 
for study with all index and reference stations marked (Figure SI 
in SI to Figure S19 in SI). 

In total, these 17 LEZs, eligible for study, contained 108 eligible 
monitoring stations with 5,3 index stations and 55 reference 
stations. The data base constructed from transferred data 
encompassed a total of 9,517,91 1 data hnes which were used as 
input to analysis. An overview is given in Table 1. 

Data analysis 

The data set structured for analysis consisted of matched 
quadruplets. A matched quadruplet comprises four pairwise 
corresponding measurement values consisting of two index- and 
two reference values. One index value and the simultaneously 
measured reference value were obtained during the active LEZ 
period, the other pair of values was obtained before introducing 
the LEZ. The pairs of values had a 364 days difference in time of 
or a multiple of 364 days, hence keeping the season, day of the 
week and time of day constant within the quadruplets. The 
allocation of reference stations to index stations was done pairwise, 
i.e., quadruplets were constructed by the data of one index station 
and allocating to it all appropriate reference stations with their 
data without a prior collapsing ("collapsing" is a technical term 
widely used in statistics describing the summary of a table in 
marginal, http://www.stata.com/manualsl3/dcollapse.pdf). The 
method has been described in detail before [45] and is a refined 
approach in comparison to other analytical strategies [46]. The 
analysis plan was critically reviewed by a chair of statistics. 

The quadruplets were analysed by the "difference score method 
in the two period case" [47]: Differences in index values were 
regressed on differences in reference values while other data were 
taken into account as covariates in fixed-effect regression analyses. 
Two types of models were fitted: a linear (additive) model and a 
log-linear (multiplicative) model. The difference of the index 
concentration data was used as the response variable in the linear 
model. The log of this response variable was entered into the log- 
linear regression model after applying an appropriate positive 
offset calculated from the data [48] . The two model types differ in 
the assumption on how covariables may influence the index station 
concentration data: on an additive scale or on a multiplicative 
scale [49,50]. 

The following covariables were taken into account in the basic 
fixed effects regression analyses: differences at reference stations in 
(ig/m^ (to control e.g. for large-scale meteorological changes and 
seasonal effects), baseline data at reference stations in \l%/'oc? (to 
control for time-dependent effects of reference data, Allison [47], 
and baseline data at index stations in |j.g/m'' (to control for 
"regression to the mean" [51]. This structure defines the basic 
regression approach. The covariables were entered into the log- 



linear (multiplicative) models after adding an appropriate offset if 
indicated [48] and then taking logs. 

The following equation describes the analysis of matched 
quadruples in the basic fixed-effect linear ("additive") regression 
model [47] 

z 

A Xnidh=E-|- ^E]j-Z|£ + bx-Xomdh,cent + bAr'Arz<ih + brTzOdh,cent + E- 
k=l 

^ ^mdh describes the difference of the index station data at 
monitoring station m between days d and d-364 ( = day d+1 in the 
year before), always at time (hour) h, i.e., xij^dh - M)Yndh (compare 
Figure 2). xomdh,cent dcuotcs the baseline value at station m on day 
d at time h, centred at the mean of all baseline values at station m. 
The terms Ar^h and r20dh,cent are the corresponding reference 
value data. The coefficient of major interest is the intercept of the 
regression model because it estimates the LEZ effect: E measures 
the mean effect across all LEZs, E+Ej^ the mean effect in zone k, 
l<k£Z. The coefficient bx accounts for "regression to the mean", 
bAr for the bias in annual levels (e.g., changed meteorological 
conditions), b^ for a time-dependent effect of reference values and 
s is the residual error of the concentration difference at the index 
stations. The second model type had the same structure but used 
logs of the terms ("log-linear", "multiplicative"). An appropriate 
small offset was added to avoid undefined logarithms [48] . 

The equation of the basic fixed-effect linear ("additive") 
regression model can be justified as follows (to keep the notation 
simple we suppress the time index: we write, eg, x^o instead of 

XzOdh). 

Let us start with assuming an ideal hypothetical situation: 
measurements are without any distortions and random errors, no 
covariates are operating. In that case we will measure the index 
concentration at time point 0 before the LEZ was introduced as a 
constant value Co at all index stations of the LEZ. After the 
introduction of the LEZ we will measure at time point 1 at the 
same index station the constant value ci. The effectiveness of the 
zone is simply E = Ci — Cq. 

But even if there are no biases, random errors and no covariates 
we do not expect to see the same E for all LEZs. The effectiveness 
may depend on characteristics of the zone k like the area of the 
LEZ, A); (e.g., we may expect a larger effect if the LEZ area is 
larger). The concentration at time point 1 may be written more 
appropriately as ci-l-fxAj^ (with a multiplicative coefficient f 
mapping the effect of the area into the concentration scale). The 
effect of zone k can be described as E = Ci+fxAj^— Co. Note that A|; 
operates as an effect modifier. 

We can take account of differences between the zones without 
referring to a specific characteristic of LEZ k, like the area: We 
may describe the effect of zone k in more abstract terms as E-l-E]^ 
[E+Ejj means E+EjjXzj^ with a multiplicative indicator Z(,, that 
takes the value 1 for zone k and zero otherwise, l<k<Z]. Ei^ ist 
the specific effect offset of LEZ k in comparison to the overall 
mean E of the LEZ effects. It is simple to extend the notation to 
cover different basefine concentrations for the different LEZs. 
Thus, E+E]j = ci+(Czi — ci)— [co+(Czo~co)] =Czi— c^o, i.e., the effect 
of zone k is the difference between the zone specific measurement 
values after (c^i) and before (Czo) the introduction of the LEZ k (ci 
and C() now denote the averages of the concentrations across all 
LEZs under study). 

Still, the approach is not very realistic. We should take into 
account background variations of the intensities, resulting from 
e.g. large area changes of the concentrations. These large area 
variations are reflected in the values r^o and r^i at the reference 
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Figure 1. Active and investigated LEZs in Germany, December 2009. The 17 LEZs included into the study are marl<ed with open red circles. 
The 1 7 LEZs that were active but were excluded according to protocol are indicated by full red circles. Capital cities of the Federal states are shown by 
black squares. This is a modification of the map as published at URL www.umweltbundesamt.de/umweltzonen. Date of access: 6th November 2009. 
doi:1 0.1 371 /journal.pone.01 02999.g001 



stations. Despite all efforts to measure the concentrations as 
precisely as possible we always will have random errors Si and Eq- 
In this extended approach the measurement values for zone k 
before introducing the LEZ are x^q = Co+(Ci,o — C())+gxr2()+S() and 
Xzi = Ci+(c2i — Ci)+gxr2i+Si after the introduction. The factor g 
measures how strong the reference values do influence the index 
values. It follows that x^i — x^o = E+E^+g* (t^i — v^q)+{&j—Zq). With 
A x^ = x^i -x^o, A r^ = r^i -r^o, s = ei "So and bAr = g we yield the 
major part of the equation of the basic frKcd-effect linear 
("additive") regression model. Note that we substituted z by m 
which means that we apply the approach in a refined way to every 



index monitoring station m. We have demonstrated above that a 
potential confounder/adjuster, like the concentration at a refer- 
ence station, enters the equation in terms of the difference of the 
values across time (e.g., A r^ = r^i — r^o). And we have seen that the 
model can be extended by potential modifiers of the LEZ effect 
("interaction terms"), like the area, by adding terms like f xA]j (in 
contrast to adjusters not as a difference in time). We wiU now 
explain in more detail why we included the effect modifying 
variables x^o and r^o additionally. 

Altman and Bland [52] and Bland and Altman [53] suggested 
including the mean value x^,,, of the concentrations at index station 
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Start of 
study 



introduction LEZ 



observation period 1 



end of 
study 



observation period 




calendar time 



XiO»in(d.1]h 



Figure 2. Index (Xjamdh) and reference concentration {r^ah) 3t index measurement station m and LZ z in observation period II with 
active LEZ (a = 1 ) and in observation period I with inactive LEZ (a = 0): matched quadruplets consisting of two index measurement 
values and two reference measurement values. The time difference between compared measurement values at day d of the year in period 2 
and d+1 of the year in period 1, always with starting time h, is not a full year but 364 days to keep the weekday constant. 
doi:1 0.1 371/journal.pone.01 02999.g002 



m as another covariate: this allows the difference A to depend 
on the average concentration at the station. This inclusion of x^,,, 
operates again a distortion due to "regression to the mean" [54]. 
This phenomenon is inevitably complicating longitudinal com- 
parisons. Baseline values that are very high due to random errors 
wiU probably not be reproduced but lower values will be 
measured, and this is so even if the nuU hypotheses of no effect 
is true [51,55-57]. A better strategy to correct for this potential 
distortion is to include x^onj^-^.n,, i.e., the baseline values at the 
index station [58,59]. Including additionally r^o, , ent was exercised 
in AUison [47], p. 10. This approach allows for a flexible 
adjustment of the annual level bias because we get rid of the 
assumption of a time-invariant effect of the reference station values 
on the index stations values. The covariates x^om, cent and r^o, cent 
are centered on the mean of the values of each measuring station 
so that the terms E und Ej; can be interpreted without further 
transformations. 

Since the impact of meteorological conditions is extremely 
relevant [e.g., 60], the following data were collected to be used in 
addition to the reference station data to control for distortions due 
to meteorological changes. We took over the height of the 
inversion base H in m, the wind velocity V in m/ s, the amount of 
precipitation P in mm/h from the PAREST project for all 
investigated measurement stations and half hours in the follow-up 
period [61-65]. 

We extended the basic regression models to the regression 
model 1 approach by adjusting additionally for the change of the 
three meteorological variables at the index stations According to 
the box model of meteorology [66] the differences were calculated 
on the additive scale after transforming the variables into 1 /H, 1 / 
(V+0.1 m/s), and 1/(P+0.1 mm/h). The smallest unit of scale was 
0.1 throughout, hence this value was used as an offset to avoid 
divisions by zero [48]. Differences were determined on the 
multiplicative scale after taking logs of these terms [48]. We 
adjusted for the time span (in years) between measurements 



considered within a quadruplet in order to adjust for trends in 
concentration levels before the LEZ was introduced. In multipli- 
cative models log of the time span was used after applying an 
appropriate offset [48]. 

In regression model 2 approach, the following time-dependent binary 
indicators were additionally adjusted for: period of school holidays (yes/ 
no), period of environmental bonus paid (yes/no) and periods when 
trucks were not allowed to enter the area where the measurement station 
was located (yes/no). In Germany a bonus was paid to car owners 
between January 14, 2009 and November 2, 2009 if they bought a new 
car with a reduced exhaust emission (littp://www.bafa.de/bafa/de/ 
wirtschaftsfoerderung/umweltpraemie/index.html). These binary indi- 
cators were entered also into the extended log-linear (multiplicative) 
models. 

This statistical approach was successfully validated in advance 
to the study in an analysis of simulated data from FU Berlin [67]. 
The simulated data was produced by the PAREST project [61,63, 
www.parest.de] . The major arm of this project was the identifi- 
cation of emission reducing strategies by simulation. Transport 
and distribution models were developed and applied, the so called 
REM-CALGRID approach [68-70]. The model was apphed to 
the city of Munich, and simulated half-hour PMiq data were 
generated for each of the five index and three reference stations 
(see Figure SIO in File SI). Data of the year 2005 were simulated 
twice, with and without adding an LEZ effect (the value of the 
imprinted effect was unknown to the analyzing working group). 
280,320 data lines were transferred. The simulated PMiq 
concentrations were analyzed and showed a mean value of about 
21 Hg/m' at the index stations and 18 |J,g/m' at the reference 
stations. The basic additive regression model estimated an LEZ 
effect of —0. 1 30 |J.g/ m' , the multiplicative model a relative change 
of —0.7%. The PAREST research report [63, www.parest.de] 
described the LEZ effect that was imprinted: PMig mean values 
are reduced in Munich city by at most 0.2 |J.g/m' or at most 1%. 
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Additive and multiplicative regression models were fitted to sub- 
sets of the data to perform sensitivity analyses: continuous 
measurement data, continuous measurement data collapsed to 

four-week averages, diffuse sampler data, pooled diffuse sampler 
and collapsed continuous data, always with and without excluding 
times with restrictions of truck traffic; quadruplets produced by 
index traffic stations only. The pooled continuous and diffuse 
sampler measurement data determined for four weeks periods was 
of major interest in this study because the annual average is the 
most critical endpoint to consider (see Introduc:tion section) and 
these data cover both types of measurement data. Because annual 
data are generally too coarse for LEZ effect estimation, we 
foUowed-up on averages over about a month. The additive and 
multiphcative regression models analysing these sets of data were 
specified with three sets of covariables as described above. All basic 
models, the models evaluating continuous data collapsed to 4-week 
averages and the models fitted to single index stations were not 
used for statistical testing of the LEZ eflFects. Tests for effects 
measured by NOg, NO, and NOx quadruplets were not 
considered as independent. According to this structure, we 
evaluated 2*3*2*2*2 = 48 statistical tests for each of the three 
endpoints. Due to this multiple testing scenario we applied an 
adapted significance level of 5%/50 = 0.1% ["family wise error 
rate", 71]. 

We fitted additionally explorative models that estimated the size 
of the LEZ effect at each index station enrolled. In addition, we 
estimated mean effects of the LEZs across the Federal states. The 
results of these exploratory analyses were mainly used for internal 
discussions of the project steering committee (see Acknowledge- 
ment). 

AU regression models used robust estimators of coefficient 
variances. AU data analyses were performed using Stata 1 1 [72] on 
a 64-bit PC. 

Results 

NO2 - continuous measurements 

The basic data consisted of 6,412,864 data lines leading to 
3,038,781 quadruplets of continuous NO2 measurement (half-hour 
averages) from 6 Federal states and 1 7 LEZs with 38 index stations 
and 45 reference stations. Table 2 gives an overview of the 
distributions observed: on average, NO2 concentrations were 
between 50 Hg/m^ and 52 Hg/m^ at the index stations and 
between 26 Hg/m^ and 27 Hg/m^ at the reference stations. The 
differences at the stations varied substantially in a range of 
hundreds of [ig/m^ upwards and downwards. A comparison of 
mean and median differences at index and reference stations 
indicated a crude LEZ effect estimate of about — 1 ng/m^. In the 
finear model 1 the absolute effect estimate was similar: —1.11 [ig/ 
(Table 3). The model 1 results showed a time-dependent 
impact of reference station data, a pronounced "'regression to the 
mean", a i k'ar influence of the three meteorological variables 
(independentiy from the crude adjustment by reference station 
data), and a downward trend of concentrations before the LEZs 
were introduced. The direction of impact of the meteorological 
variables was as expected: the smaller H, V, or P die larger the 
index N02 concentrations. In linear model 2 the LEZ effect 
estimate was slightly more pronounced: — 1.85 |J,g/m3. In the log- 
finear model 1 (multiplicative approach), the relative effect 
estimate was 0.979, i.e., a reduction of 2.1% was found (Table 4). 
The estimated impact of covariables agreed with the finding in the 
corresponding linear model. When applying regression model 2 
the relative LEZ effect estimate was 0.961, i.e., the reduction was 
estimated to be 3.9%. 



NO2 - pooled continuous and diffuse sampler 

measurement data 

6,133 data lines and 4,095 quadruplets of NO2 pooled 
continuous and diffuse sampler measurement data (averaging 
period: four weeks) were examined from 17 LEZs with 53 index 
stations and 55 reference stations. A crude comparison based on 
the observed distributions revealed a LEZ effect of about — 
0.2 |tg/m" to -0.6 ng/m' (Table 5). U sing the linear model 1 
approach the absolute effect estimate was —0.826 pg/m'^ 
(Table 6). The meteorological variables showed no substantial 
impact due to the long averaging period. Model 2 estimated the 
LEZ effect as —1.73 |ig/m'^. The log-linear modelling led to a 
relative effect of 0.980 (Table 7, model 1) or 0.961 (model 2). 
Table S 1 in File S 1 provides a detailed overview of the results 
when fitting a series of models to analyse the NO2 measurements. 
LEZ effect estimates were about —1 |tg/m^ to — 2 |J.g/m^ 
(additive models) or —2% to —4% (multiplicative models). 

NO - pooled continuous and diffuse sampler 
measurement data 

A total of 5,790 data lines from 17 LEZs with 46 index stations 
and 54 reference stations were available to analyse pooled 
continuous and diffuse sampler NO measurement data. A 

descriptive analysis of the 4,005 quadruplets indicated a LEZ 
effect of about 0 |tg/m'^ to —1 pg/m"* (Table 8). Using the 
additive approach the absolute effect estimate was —1.13 pg/ m ^ 
in model 1 (Table 9). When the model specification 2 was appfied 
the LEZ effect estimate changed the sign: -1-0.38 |a,g/m^, i.e., no 
reduction was indicated in this extended model type. The log- 
linear regression model of type 1 yielded a relative effect estimate 
of 0.968 (Table S2 in File SI). The direction of the estimated 
relative effect changed when model 2 was applied: +1.20. 

NOx - pooled continuous and diffuse sampler 
measurement data 

The analysis of pooled continuous and diffuse sampler NO, 
measurement data was performed using 4,005 quadruplets that 
originated from a set of 5,790 data lines generated by 46 index 
stations and 54 reference stations of 17 LEZs. According to the 
distributions of differences a crudely estimated LEZ effect (based 
on averages or medians) was present of about —0.2 pg/m'' to — 
1.3 \lg/m^ (Table S3 in File SI). Adjusting for covariables in linear 
model 1 returned an absolute effect estimate of —1.74 Hg/m^ 
(Table S4 in File SI). The adjustment for further covariables 
(regression model 2) led to an effect estimate of —0.89 pg/m'^. 
When the k)g-linear model 1 was used (Table S5 in File SI), a 
relative LEZ effect of 0.976 was found. The adjustment for 
additional covariates (model 2) led to a change in direction: the 
relative effect was estimated as 1.048. 

Summary of Results for NO2, NO, and NOx 

Table 10 gives an overview of the findings for NO2, NO, and 
NOx. The mean concentration levels at the index stations were 
about 50 Hg/m^ for NO2 and for NO, thus, about 100 Jig/m^ for 
NOx. Model 1 analyses showed reductions of the concentrations 

after introducing the LEZs. Although small, all effect estimates 
were statistically significant at the 0. 1 % level. Model 1 estimates 
based on an additive structure gave compatible findings to the log- 
Unear multiplicative approach (e.g., 2% of 50 (tg/m'*= 1 (tg/m'*). 
The model 1 LEZ effect estimates were similar to, but slightiy 
more pronounced than crude LEZ effect estimates based on direct 
comparisons of the measurement differences at index stations and 
reference stations within the quadruplets while ignoring the impact 
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Table 2. NOj: Quadruplets of continuous N02-measurements: index stations (Ind), reference stations (Ref) before (pre) and after 
(post) introduction of LEZ. 





Statistic 


lnd,pre 


lnd,post 


Ref,pre 


Ref.post 


ind.dlff. 


Ref.dlff 


N 


3038781 


3038781 


3038781 


3038781 


3038781 


3038781 


min 


0.4 


1.3 


0.4 


0.5 


-330 


-215 


p5 


12 


11 


4.0 


4.0 


-52 


-33 


p50 


45.0 


43.8 


20.7 


20.6 


-1.0 


0.0 


mean 


51.959 


50.831 


26.383 


26.17 


-1.128 


-0.212 


p95 


115 


114 


68 


67 


49 


33 


max 


392 


436 


248 


434 


375 


317 



Ind.diff and Ref.diff denote differences between index measurements and between reference measurements (negative post-pre differences indicate lower values after 
introduction of LEZ). Concentrations measured in fig/m^. 

N: number of quadruplets; Min: minimum, p5: 5**^ percentile, p50: median, mean: arithmetic average, p95: 95*"^ percentile, max: maximum, unit: 
doi:1 0.1 371/journal.pone.01 02999.t002 



of covariables (compare Tables 2, 5, and 8 and Table SI in SI). 
AU analyses point to the conclusion that on average the 
concentration reducing elTect of LEZs was smaller than 2 |J.g/ m 
for each of the three components NO2, NO, and NO^, i.e., not 
higher than about 4%, when considering all investigated index 
stations. However, breaking down the analyses by Federal states or 
LEZs yielded heterogeneous estimates of effects. 

The NO2 analysis was based on 192 comparisons of index vs 
reference stations, among them were 3 1 index stations character- 
ized as "background", one characterized as "industry" and 160 as 
"traffic" stations. We performed a sensitivity analysis by restricting 
the evaluation to the stations close to traffic. The additive linear 
type 2 model estimated an effect of —1.73 |J.g/ m at all index 
stations (see last line in Table SI in SI). When the analysis only 
accounted for the traffic stations we got a slightly more 
pronounced LEZ effect estimate of —2.26 |ig/m'' (3,406 quadru- 
plets, pooled data: four week averages). An analysis of the 
continuous data yielded almost the same result: —2.35 M-g/m'' 
(2,105,702 quadruplets, half-hour averages). 

Discussion 

In this study we analysed the effect of introducing the "Low 
Emission Zone (LEZ) of pollutant group 1" (which restricts from 



entering Diesel cars of an European emission standard below Euro 
2 without particulate reduction) on NO2, NO, and NOx 
concentrations in Germany. We included as many LEZs as 
possible (17 out of 34 in 2009 met our inclusion criteria) into a 
homogeneous analysis of nitrogen oxide data measured before and 
after the introduction of LEZs of pollutant group 1 until the end of 
2009. We used matched quadruplets of index and reference station 
values and analysed the changes in concentrations with frKcd-effect 
regression models while adjusting for important covariables. We 
performed sensitivity analyses by applying two model structures 
(additive and multiplicative) with varying sets of covariables to 
different subsets of the data. We based our study on precisely 
matched quadruplets to avoid distortions and to increase validity. 
A potential downside of the increased validity is a loss in precision 
due to the reduced data set eligible for analysis. However, the loss 
in power was negligible in this application because P-values were 
small even when taking multiple testing into account [73]. The 
statistical approach was successfully validated in advance to the 
study in an analysis of simulated data from FU Berlin [67]. We 
checked whether the adjustment in one model that analyzed all 
LEZs simultaneously and assumed unknown but identical 
covariate coefficients was appropriate for all LEZs. To do so we 
evaluated each LEZ separately and performed a meta-analysis on 
the findings. The precision weighted mean of the effect estimates 



Table 3. NO2: Linear (additive) model 1 evaluating the quadruplets of continuous N02-nneasurennents. 





Ind.diff 


Coef. 


Std. Err. 


t 


P 


95% Conf. 


Interval 


Refdiff 


0.677 


0.001 


623 


<0.001 


0.675 


0.679 


Ref base 


0.509 


0.001 


417 


<0.001 


0.507 


0.512 


Ind. base 


-0.644 


0.001 


-845 


<0.001 


-0.645 


-0.642 


Diff 1/H 


564 


2.8 


200 


<0.001 


558 


569 


Diff 1/V 


3.12 


0.025 


126 


<0.001 


3.07 


3.17 


Diff 1/P 


0.088 


0.004 


23.6 


<0.001 


0.081 


0.095 


Time.diff 


-0.399 


0.023 


-17.7 


<0.001 


-0.443 


-0.354 


E 


-1.112 


0.013 


-87.3 


<0.001 


-1.137 


- 1 .087 



Regression coefficient, robust standard errors of coefficient, t-statistic, two-sided P-value, and 95%-confidence interval of coefficient. The absolute LEZ effect estimate is 
given by the coefficient E in fig/m^ (<0: concentration is lowered by LEZ). 

Covariables: difference in reference stations Ref.diff in |Ig/m^ centered reference baseline concentration Ref.base in ^g/m^ centered index baseline concentration 
Ind. base in p.g/m^ Diff 1/H - difference in 1 /(height of inversion layer) in 1/m, Diff 1/V = 1/{wind velocity+O.OI) in (m/s)~\ difference in 1/P = 1 /(amount of precipitation-l- 
0.01) (mm/h)~\ centered difference in time Time.diff in years. 
doi:l 0.1 371/journal.pone.01 02999.t003 
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Table 4. NO2: Log-linear (multiplicative) model 1 evaluating the quadruplets of continuous NOj-measurements. 





In Ind.diff 


Coef. 


Std. Err. 


t 


P 


95% Conf. 


Interval 


In Refdiff 


0.296 


0.001 


582 


<0.001 


0.295 


0.297 


In Ref.base 


0.265 


0.002 


141 


<0.001 


0.262 


0.269 


In Ind.base 


-0.981 


0.002 


-412 


<0.001 


-0.986 


-0.977 


Diff In 1/H 


0.071 


0.001 


194 


<0.001 


0.070 


0.071 


Diff In 1/V 


0.12 


0.001 


301 


<0.001 


0.119 


0.121 


Diff In 1/P 


0.014 


0.001 


29 


<0.001 


0.013 


0.015 


In TIme.diff 


-0.025 


0.001 


-32 


<0.001 


-0.027 


-0.024 


In E 


-0.021 


0.0 


-77 


<0.001 


-0.022 


-0.021 



Regression coefficient, robust standard errors of coefficient, t-statistic, two-sided P-value, and 95%-confidence interval of coefficient. The relative LEZ effect estimate is 
given by the coefficient E (<1: concentration is lowered by LEZ). 
relative effect E = 0.979, 95% Conf Interval = 0.979, 0.980. 

Covariables, before taking logs: difference in reference stations Ref diff in fig/m', centered reference baseline concentration Ref base in )ig/m^, centered index baseline 
concentration Ind.base in jig/m^, 1/H = 1 /(height of inversion layer) in 1/m, 1/V = l/(wind velocity+0.01) in (m/s)"\ 1/P = 1 /(amount of precipitation+0.01) (mm/h)"\ 
centered difference in time TIme.diff in years. 
doi:l 0.1 371 /journal.pone.Ol 02999.t004 



at all iiidex stations (n= 192, —1.71 [ig/m ) was almost identical 
to the overall additive linear type 2 model (—1.73 |J.g/m'', see last 
line in Table SI in SI). We conclude that the fitted single model 
that evaluated all LEZs simultaneously was appropriate and did 
not suffer fi-om an insufficient adjustment. 

As an overall finding the average effect of LEZ introduction on 
nitrogen oxide concentrations (NO2, NO, and NO^ = NO2+NO) 
was not higher than 2 |a,g/m^ at all index stations, i.e., not higher 
than about 4%. The effect was only slightly larger when we 
restricted the analyses to stations close to traffic. In the main 
analyses the coefficients describing the reductions were statistically 
significant on the 0. 1 % level, i.e., after taking multiple testing into 
account. We note, however, that the P-values calculated are 
potentially too small because autocorrelations in the data were not 
taken into account. 

We detected a substantial heterogeneity of effects across the 
investigated LEZs and Federal states. However, this finding is not 
surprising because 

— the realisation of LEZs differed between states and within states 
(e.g., date of introduction, covered population and area of 
LEZs differ (compare Table 1), some operate together with an 
additional restriction of van traffic) 



- the degree of representativeness of monitoring stations inside 
the LEZs differs across LEZs (index stations: distances from 
centre/border of LEZ differ, used as background or hot spot 
stations and sometimes placed in street canyons) 

- the degree of representativeness of monitoring stations outside 
the LEZs differs across LEZs (reference stations: distances from 
LEZ differ, traffic conditions differ) 

- the applied measuring systems differ (continuous chemUumi- 
nescense procedure vs diffuse long-term sampling with 
chromatography). 

The large variation of LEZ effect estimates across the LEZs 
should be put into perspective by considering the phenomenon of 
"regression-to-the-mean" [51]. Due to this phenomenon we 
expect that single observations with high baseline values show 
potentially decreasing trends - and low baseline values potentially 
increasing trends. This is true even under the null hypothesis of no 
causal LEZ effects on nitrogen oxide concentrations. "Regression- 
to-the-mean" has been shown to be rather pronounced in this 
study. Thus, the interpretation of single LEZs effect estimates is 
clearly limited and we will not report any details with the consent 
of the involved state institutions who performed the measurements 
(see Acknowledgements). 



Table 5. NO2: Quadruplets of pooled continuous and diffuse sampler N02-measurements: index stations (Ind), reference stations 
(Ref) before (pre) and after (post) introduction of LEZ. 



Statistic 


Ind.pre 


Ind.post 


Ref.pre 


Ref.post 


Ind.diff. 


Ref.diff 


N 


4095 


4095 


4095 


4095 


4095 


4095 


min 


13 


11 


6 


5 


-32 


-23 


P5 


26 


24 


n 


10 


-14 


-13 


p50 


57.31 


54.23 


41.04 


39.78 


-2.34 


-1.71 


mean 


56.298 


54.246 


40.007 


38.18 


-2.053 


- 1 .824 


p95 


84 


82 


66 


64 


12 


8 


max 


134 


136 


75 


76 


41 


25 



Ind.diff and Ref diff denote differences between index measurements and between reference measurements {negative post-pre differences Indicate lower values after 
Introduction of LEZ). Concentrations measured in fig/m^. 

N: number of quadruplets, MIn: minimum, p5: 5**^ percentile, p50: median, mean: arithmetic average, p95: 95**^ percentile, max: maximum. 
doi:l 0.1 371/journal.pone.Ol 02999.t005 
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Table 6. NOy. Linear (additive) model 1 evaluating the quadruplets of pooled continuous and diffuse sampler NOj-measurements. 



Ind.diff 


Coef. 


Std. Err. 


t 


P 


95% Conf. 


Interval 


Ref.diff 


0.564 


0.018 


31.7 


<0.001 


0.529 


0.599 


Ref.base 


0.394 


0.025 


15.9 


<0.001 


0.345 


0.442 


Ind.base 


-0.695 


0.019 


-36.7 


<0.001 


-0.732 


-0.657 


Diff 1/H 


-65.1 


20.8 


-3.13 


0.002 


-106 


-24.3 


Diff 1/V 


0.034 


0.148 


0.23 


0.816 


-0.256 


0.325 


Diff 1/P 


0.009 


0.028 


0.30 


0.764 


-0.047 


0.064 


Time.diff 


0.369 


0.189 


1.95 


0.051 


-0.001 


0.739 


E 


-0.826 


0.123 


-6.71 


<0.001 


- 1 .068 


-0.585 



Regression coefficient, robust standard errors of coefficient, t-statistic, two-sided P-value, and 95%-confidence interval of coefficient. The absolute LEZ effect estimate is 
given by the coefficient E In fig/m^ (<0: concentration Is lowered by LEZ). 

Covariables: difference In reference stations Ref.diff In |ig/m^, centered reference baseline concentration Ref base in ^g/m^, centered Index baseline concentration 
Ind.base in ^g/m^, Diff 1/H = difference in 1 /(height of Inversion layer) in 1/m, Diff 1/V = 1/{wind veloclty+0.01) in (m/s)"\ difference in 1/P = 1 /(amount of preclpitation+ 
0.01) (mm/h)"\ centered difference in time Time.diff in years. 
doi:l 0.1 371/journal.pone.Ol 02999.t006 



Models of type 2 showed more instability and returned positive 
effect estimates in some situations (see Table 10). Regression 
model 2 included as additional variables time-dependent binary 
indicators for period of school holidays, period of environmental 
bonus paid and periods when trucks were not allowed to enter the 
area where the measurement station was located. In some LEZs 
these variables were highly correlated with the active LEZ periods 
so that unstable findings due to coUinearities can be expected. 
Such coUinearities can introduce a bias away from the nuU and 
may generate exaggerated negative or positive model coefficients 
even if the true effects are near to zero [74]. Log-linear models 
showed to be more sensitive to these distortions. This may indicate 
a less appropriate modelling of the data when assuming 
multiplicative effects of covariates. 

There are evaluations available concerning potential effects of 
LEZs on NO2 concentrations summarized by the German Federal 
Environmental Agency [7]: A total NO2 reduction by 5% and a 
local traffic-related NO2 reduction by 1 2 % may be reached given 
"LEZ of pollutant group 3" so that only cars with a green sticker 
(Diesel vehicles of Euro 6, 5, 4 or Euro 3 with particle fdter, 



gasoline cars with catalytic converter) are allowed to enter the LEZ 
[3] . This statement is based mainly on preliminary evaluations of 
the Berlin LEZ data by Rauterberg-Wulff and Lutz [29]. Puis and 
Jager-Ambrozewicz [75] reported for the Frankfurt LEZ and an 
observation period until the end of 2011 effects of less than 3% 
which is closer to our present findings although they also cover a 
period of "LEZ of pollutant group 2" after Janl, 2010. Only cars 
with a yellow sticker (Diesel vehicles of Euro 3 or 4 standard or 
Euro 2 with particle filter, gasoline cars with catalytic converter) 
were allowed to enter the Frankfurt LEZ after Jan 1 , 2010 [3] . 
Bruckmann et al. [6] reported reductions of the annual average of 
NO2 concentrations up to 2% associated with the introduction of 
"LEZs of pollutant group 1" in North-Rhine Westphalia, and an 
absolute LEZ effect of about — 1 .2 \lg/ ni' . In Hannover no NO2 
reduction could be shown after introducing an "LEZ of poUutant 
group 1" [76]. All of these statements, however, were based on 
crude comparisons without suflTiciendy adjusting for important 
covariates like weather conditions, and traffic restrictions etc. Only 
Puis and Jager-Ambrozewicz [75] applied a more sophisticated 
approach. They performed a time-series analysis and fitted 



Table 7. NO2: Log-linear (multiplicative) model 1 evaluating the quadruplets of pooled continuous and diffuse sampler NOj- 
measurements. 



in Ind.diff 


Coef. 


Std. Err. 


t 


P 


95% Conf. 


Interval 


In Refdiff 


0.359 


0.015 


23.9 


<0.001 


0.329 


0.388 


In Ref.base 


0.024 


0.005 


4.49 


<0.001 


0.014 


0.034 


In Ind.base 


-0.039 


0.004 


-10.1 


<0.001 


-0.047 


-0.032 


Diff In 1/H 


-0.016 


0.003 


-5.32 


<0.001 


-0.022 


-0.01 


Diff In 1/V 


0.013 


0.003 


4.21 


<0.001 


0.007 


0.02 


Diff In 1/P 


0.011 


0.003 


3.68 


<0.001 


0.005 


0.017 


In Time.diff 


0.005 


0.006 


0.95 


0.34 


-0.006 


0.017 


In E 


-0.02 


0.003 


-7.82 


<0.001 


-0.025 


-0.015 



Regression coefficient, robust standard errors of coefficient, t-statistIc, two-sided P-value, and 95%-confidence interval of coefficient. The relative LEZ effect estimate is 
given by the coefficient E (<1: concentration is lowered by LEZ). 
relative effect E = 0.980, 95% Conf Interval = 0.975, 0.985. 

Covariables, before taking logs: difference in reference stations Refdiff in ^g/m^, centered reference baseline concentration Ref base in \ig/m^, centered Index baseline 
concentration Ind.base In jig/m^, 1/H = 1 /(height of Inversion layer) in 1/m, 1/V = l/(wind velocity-t-0.01) in (m/s)"\ 1/P = 1 /(amount of precipltatlon+0.01) (mm/h)"\ 
centered difference in time Time.diff In years. 
doi:l 0.1 371/journal.pone.Ol 02999.t007 
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Table 8. NO: Quadruplets of pooled continuous and diffuse sampler NO-measurements: Index stations (Ind), reference stations 
(Ref) before (pre) and after (post) introduction of LEZ. 





Statistic 


Ind.pre 


Ind.post 


Ref.pre 


Ref.post 


Ind.diff. 


Ref.diff 


N 


4005 


4005 


4005 


4005 


4005 


4005 


min 


2.7 


2.5 


1.0 


1.0 


-86 


-50 


p5 


4.6 


4.4 


2.2 


2.2 


-22 


-21 


p50 


48.54 


43.44 


31.83 


26.32 


-2.59 


-1.26 


mean 


49.479 


46.373 


34.153 


31.025 


-3.105 


-3.129 


p95 


95 


95 


86 


81 


17 


13 


max 


230 


251 


139 


104 


77 


56 



Ind.diff and Ref.diff denote differences between index measurements and between reference measurements (negative post-pre differences indicate lower values after 
introduction of LEZ). Concentrations measured in fig/m^. 

N: number of quadruplets, Min: minimum, p5: 5**^ percentile, p50: median, mean: arithmetic average, p95: 95**^ percentile, max: maximum. 
doi:l 0.1 371/journal.pone.Ol 02999.t008 



regression models for the Frankfurt LEZ. These models', however, 
were not correctly specified as they did not include difTerences of 
the covariables but the absolute values only, and so they could not 
control for potential confounding effects although this was 
intended by the authors. AH publications cited above reported 
only on individual LEZs or certain Federal states in Germany and 
not on the LEZ effect on the national level. Generalisations from 
these data are problematic because of the heterogeneous 
configurations of LEZs. A realistic estimate should be based on 
a homogeneous analytical approach covering as many LEZs and 
Federal states simultaneously as possible, as performed in this 
study. 

Table 1 1 presents an overview of other study results published 
in the peer-reviewed literature on forecasted or measured LEZ 
effects on NO2 concentrations. 

Our results are in good accordance with the prognosis study 
PAREST of FU Berlin [61,63]. An extensive description of the 
project is available [77]. The prognoses of PAREST are 
comparable with our estimates at all index stations because 
PAREST worked with an area coarseness defined by grid square 
of about 1 km x 1 km and, thus, cannot estimate changes at single 
stations. Duyzer et al. [78] studied whether monitoring station 
data are representative for the population living in the area and 
concluded that the background station data are more appropriate 



to describe the impact on the citizens than the hot spot traffic 
stations. We conclude that the findings of PAREST and our results 
about the effect at all index stations should be preferred in an 
evaluation (not the effect estimates restricted to the traffic stations). 
PAREST predicted LEZ effects on NO2 levels assuming that only 
cars with green stickers are allowed to enter (LEZ of pollutant level 
3). For the Berlin LEZ the authors calculated a reduction of about 
1 |Xg/m' to 1.3 Hg/m"* in the city centre (relative: 3% to 5%), for 
the Munich LEZ a reduction of 1 \xg/m'^ in the city centre 
(relative: up to 5%), for the Ruhr area a reduction of 1 |Xg/m'^ to 
1.7 (ig/m'^ (relative: 3% to 4%). Setting the whole Ruhr area to a 
LEZ of pollutant level 3 lead to the prognosis of a reduction in 
NO2 concentrations of 1 Hg/m'^ to 2 Hg/m'' (relative: 3% to 6%). 
It needs to be taken into account that these prognoses by PAREST 
are based on the pollutant level 3 LEZ scenario. We do not expect, 
therefore, that our findings from this study may change relevantiy 
if the LEZs are extended to cover larger areas or if stricter traffic 
restrictions are applied. 

A very large LEZ was introduced in London as a congestion 
charging zone. However, only prognoses of the potential LEZ 
effect on nitrogen oxide concentrations are available. NOx 
reductions between 3.8% in 2008 up to 7.3% in 2012 along 
roadways were predicted in a modelling scenario for the London 
LEZ with vehicles and buses required to meet Euro 4 standards 



Table 9. NO: Linear (additive) model 1 evaluating the quadruplets of pooled continuous and diffuse sampler NO-measurements. 





Ind.diff 


Coef. 


Std. Err. 


t 


P 


95% Conf. 


Interval 


Refdiff 


0.666 


0.022 


30.5 


<0.001 


0.623 


0.708 


Ref.base 


0.496 


0.026 


19.1 


<0.001 


0.445 


0.547 


Ind. base 


-0.473 


0.022 


-21.5 


<0.001 


-0.516 


-0.430 


Diff 1/H 


75.8 


29.5 


2.57 


0.010 


18.0 


134 


Diff 1/V 


-0.571 


0.216 


-2.65 


0.008 


-0.994 


-0.148 


Diff 1/P 


0.028 


0.048 


0.59 


0557 


-0.066 


0.123 


Time.diff 


0.245 


0.313 


0.78 


0.433 


-0.369 


0.860 


E 


-1.128 


0.218 


-5.19 


<0.001 


-1.555 


-0.702 



Regression coefficient, robust standard errors of coefficient, t-statistic, two-sided P-value, and 95%-confidence interval of coefficient. The absolute LEZ effect estimate is 
given by the coefficient E in fig/m^ (<0: concentration is lowered by LEZ). 

Covariables: difference in reference stations Ref.diff In |Ig/m^ centered reference baseline concentration Ref.base in ^g/m^ centered index baseline concentration 
Ind. base in p.g/m^ Diff 1/H - difference in 1 /(height of Inversion layer) in 1/m, Diff 1/V = 1/{wind veloclty+O.OI) in (m/s)~\ difference in 1/P = 1 /(amount of preclpitation-l- 
0.01) (mm/h)~\ centered difference in time Time.diff in years. 
doi:l 0.1 371/journal.pone.01 02999.t009 
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Table 10. Summarized Results on Nitrogen Oxides. 











IVIodei 1 Estimates 




IVIodei 2 Estimates 






Evaluation 
Period 


Mean index Concentration/|ig/m^ 


Additive Effect/^g/m^ 


IVIultipllcative Effect/% 


Additive Effect/|ig/m^ 


Multiplicative 
Estimate/% 


N02 


0.5 h 


50 


-1.1 


-2.1 


-1.9 


-3.9 




4 weeks 


55 


-0.8 


-1.9 


-1.7 


-3.9 


NO 


4 weeks 


48 


-1.1 


-3.2 


+0.4 


+1.9 




4 weeks 


103 


-1.7 


-2.4 


-0.9 


+4.8 



Model 1 covariables: difference in reference stations Refdiff in ^g/m^ centered reference baseline concentration Ref.base in ^ig/m^, centered index baseline 
concentration Ind.base in [ig/rvi^, Diff 1/H = difference in 1 /(height of inversion layer) in 1/m, Diff 1/V= 1/(wind veloclty+0.01} in (m/s)"\ difference in 1/P = 1 /(amount of 
precipitation+O.Ol} (mm/h)"\ centered difference In time Time.diff in years. 
doi:l 0.1 371 /journal.pone.Ol a2999.tai 0 



compared to current LEZ restrictions for Euro 3 vehicles [46]. prognoses published by Tonne et al. [79], who estimated for the 

The authors stated that despite of the large area of the London London Congesting Charge Scheme a small decrease for NO2 of 

LEZ, the predicted changes in NO2 (and PMio) were generally —0.64 |J.g/m^ only, corresponding to —1.1%. 
small. Their modelled results stay partly in contrast to the 



Table 11. Overview of studies estimating the effect of LEZs on NO2 concentrations. 



Study 


Country 


City/Area 


Intervention 


NO2 effect/fig/m^ 


Comments 


Builtjes et al. 2012, 


Germany 


Berlin 


LEZ of level 3: forbidden<Euro 4 


-1.3 




Stern 2013 




Munich 




-1.0 








Ruhr Area 




-1.7 




Tonne et al. 2008, 


UK 


London 


Congestion Charging Zone: 
forbidden<Euro 4 


-0.64 (-1.1%) 


simulation study 2008 


Kelly et al 2011 








reduction up to —7.3% 


simulation study 201 1 


Briggs 2008 


Italy 


Rome 


2 LEZs: 


-2.3/-3.0 


simulation study: main effect 
due to the exclusion of Euro 0 
vehicles 


Euro 0 forbidden 2 LEZs: 


Cesaroni et al. 
2012 


Italy 


Rome 


forbidden<Euro 4 


-3.0/-4.1 




Boogaard et al. 
2012 


The Netherlands 


Amsterdam 


LEZ: 


-4.5 


analysis of measurements: 
crude comparisons, no 
covariates taken into account 






Den Bosch 


trucks forbidden 


(not statistically significant} 








The Hague 


<Euro 2 and Euro 3 trucks only 
allowed if retrofitted 






Tilburg 


Utrecht 


Johansson et al. 
2009 


Sweden 


Stockholm 


Road Pricing System: Vehicles 
travelling into and out of the 
charge cordon were charged for 
every passage during weekdays 


NO,: 


analysis of measurements: 
crude comparisons, no 
covariates taken into account 


-0.23 (Greater Stockholm) 


-0.81 (inner city) 


This study 


Germany 


17 cities 


LEZ of level 1: Diesel passenger 
cars, trucks and buses 
forbidden<Euro 2 without 
particulate reduction system 


NO2, NO, NO,: 


analysis of measurements: 
construction of matched 
quadruplets, regression 
analyses with covariates, all 
estimated effects statistically 
significant at the 0.1% level 



reduction less than —2.0 
(-4%) 

[PMio: reduction less than - 
0.2 |ig/m^ (-1%)] 
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The INTARESE project [80] modeled NO^ concentration 
changes for both LEZs in Rome and confirmed this finding of only 
small additional gains by stricter traffic restrictions. The main 
reductions were expected to be achieved already by excluding 
Euro 0 cars: —2.3 |ig/m^ or —3.0 (tg/m\ If only Euro 4 car.s were 
allowed to enter the LEZs the reductions were expected to 
increase only slightly to -3.0 |tg/m^ or —4.1 \lg/m^ [81]. 

The "Stockholm Trial" involved a road pricing system to 
improve the air quality and reduce traffic congestion. The test 
period of th(- trial was Januar)- 3, 2006 to July 31, 2006. Vehicles 
travelling into and out of the charge cordon were charged for 
every passage during weekdays. Annual mean contributions to 
total levels of nitrogen oxides from emissions from road traffic with 
and without charges according to the Stockholm Trial were 
estimated. NO^ concentrations were lowered in periods with 
charges, but the study showed a small decrease only: —0.23 |Xg/ 
m'^ (Greater Stockholm) and -0.81 Hg/m'^ (inner city) [82]. No 
multivariable modeling was tried. 

Boogaard et al. [83] analyzed measurements of NO2 and NO^ 
conducted simultaneously at eight streets, six urban background 
locations and four suburban background locations before (2008) 
and two years after implementation of an LEZ (20 1 0) in five cities 
of The Netherlands (8 index stations, 4 reference stations). Index 
concentrations were lower in 2010 than in 2008 (NO2: —4.5 ng/ 
m^, NO^: —6.1 Hg/m^) but the differences were not statistically 
different. The study performed only crude comparisons and did 
not apply regression techniques to adjust for covariables. 

The present study can be regarded as one of the most 
comprehensive approaches so far, analysing measurement data 
of nitrogen oxides concentrations in order to assess LEZ effects. 
The LEZ pollutant group 1 reduction eflFect on nitrogen oxides 
(NO2, NO, and NOJ was estimated as being no higher than 2 |a.g/ 
m' at all index stations and index traffic stations, i.e., no higher 
than about 4'/o. This estimate based on measurement data can be 
rated as the most profound currently available. This result also 
needs to be interpreted in the light of the existing EU limit values 
because LEZs are often supposed to be the most effective measure 
that cities can take to reduce air pollution problems in their area 
[84]. The respective NO2 concentration limit [24] enforced in 
Germany since 2010 is 40 |tg/m'^ (1 year average). Values are in 
excess and about 69% of all German traffic stations showed 
annual averages higher than 40 Hg/m'^ [2.5]. The four week 
averages of NO2 concentrations at index stations after LEZ 
introduction were found to be 55 |J.g/m' (median and mean) or 
82 |tg/m'^ (95th percentile). It follows that the estimated reduction 
of NO2 concentrations in the range of 2 Hg/m'' appears to be of 
negligible impact when the current concentration levels should be 
lowered to the EU limit. The same judgement seems to apply on 
the EU level where the NO2 concentrations were reported to show 
a pronounced excess in many cities [26]. 

Regarding the information from the HBEFA [85] for real 
driving conditions in Germany, Austria and Switz(;rland with 
respect to vehicles that meet Euro 5 and 6 emission standards, no 
noteworthy reductions of NO2 and NO^ immissions are to be 
expected until a remarkable share of vehicles with NOx after 
treatment systems (Euro 5 for HD trucks and Euro 6 for passenger 
cars) will be on the street [85]. 

The Handbook of Emission Factors for Road Transport 
(HBEFA) was originally developed on behalf of the Environmental 
Protection Agencies of Germany, Switzerland and Austria. In the 
meantime, further countries (Sweden, Norway, France) as well as 
the JRC (European Research Center of the European Commis- 
sion) are supporting HBEFA. HBEFA provides emission factors, 
i.e. the specific emission in g/km for all current vehicle categories 



(PC, LDV, HDV, buses and motor cycles), each divided into 
different categories, for a wide variety of traffic situations (http:// 
www.hbefa.net/ e /index, html). 

Interestingly, remarkable differences in NO^ and NO2 emissions 
from passenger cars and light duty vehicles are documented when 
low test cycle emissions were compared with relatively higher 
N0x/N02 concentrations measured along roadsides [86,87]. 

We analysed PMio concentrations additionally [32] from 19 
German LEZs. From about 2005 until the end of 2009 continuous 
half-hour measurement values as well as gravimetrically deter- 
mined daily measurements of PMio were collected. Two 
continuous procedures were used to measure mean PMio 
concentrations per half-hour intervals [38,88]: 

- Absorption of P-radiation (BA). The particulate matter is 
deposited on a filter tape and the change in P-ray transmission 
is measured. 

- Tapered Element Oscillating Microbalance (TEOM). An 
inertial balance directly measures the mass collected on an 
exchangeable filter cartridge by monitoring the corresponding 
frequency changes of a tapered element. 

In addition, gravimetric samplers were used to measure daily 
averages of PMk, concentrations [49,88,89]. 2,110,803 quadru- 
plets of continuous PMio and 15,735 gravimetric quadruples were 
identified leading to 61,169 quadruplets based on daily PMio 
averages. The analyses showed that best LEZ effect estimates were 
^0.2 |J.g/m'^ at all index stations, i.e., the relative PMio reduction 
^ 1 % . Best estimates at all index stations near traffic (excluding 
urban background and industry index stations) were below 1 \lg/ 
m"* (less than 5%, resp). Effects were smaller than predicted prior 
to the introduction of LEZs. Limited data (1750 quadruplets of 
monthly averages) were also available to estimate the effects on 
soot parameters (elemental carbon, organic carbon and total 
carbon). The average of total carbon concentrations was estimated 
as 13 |ig/m' and LEZ effect estimates were about —0.55 \lg/in 
or —4.2%. For PM2.5 only 650 quadruplets based on half-hour 
data and 99 quadruplets of daily concentration averages could be 
analyzed. The PM2.3 concentration mean was found at 17 |a,g/m^. 
All LEZ effect estimates on PM2 were positive, i.e., no indication 
of reduced concentrations after the introduction of the LEZs was 
found. 

Due to the proven marginal reduction of nitrogen oxide 
concentrations (NO, NO2, NO^), LEZ as a regulatory action 
cannot be seen as an efficient measure to substantially reduce 
ambient nitrogen oxide exposures in the cities. Beyond that, this 
result is in good accordance to the effectiveness of LEZs on the 
reduction of PMio, too [32]. As predicted [33], long-term 
compliance problems with ambient air NO2 concentrations should 
be expected even if LEZs were introduced or enlarged for the 
purpose of NO2 reductions in cities. 

The approach can be extended to account for other variables 
that are considered relevant [45]. Such data can only be used if 
these data are homogeneously available at all index and reference 
stations and are also available before and after the introduction of 
LEZs. Traffic density and car fleet properties are such variables of 
interest that do not meet the inclusion criteria: there are almost no 
data available in Germany to describe differences in flow of traffic 
and car fleet properties between index and reference stations and 
across time. To put this into perspective, we like to note first that 
changes in traffic density and car fleet properties are potentially 
affected by LEZs. It follows that traffic density and fleet properties 
should be considered as potential outcomes of LEZ introduction 
and not only as confounders of LEZ effects. This means that these 
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data must not be accounted for by covariables in regression 
modelling even if the data were available in such a way that the 
inclusion criteria were met. Anyhow, authors who described 
changes in traffic-flow in Berlin argued against the interpretation 
that LEZs caused such displacements of traffic-flow from inside the 
LEZ to the reference stations [31]. S(;cond, we note that the 
missing information on traffic density and fleet properties can be 
used to argue for biases in both directions. On the ones side, traffic 
could be displaced from the LEZ area to the reference stations 
outside so that the concentrations are underestimated inside but 
overestimated outside the LEZ, causing a potential overestimate of 
the LEZ effect. On the other side, if the car fleet is renewed not 
only inside the LEZ but also outside at the reference stations this 
may lead to a potential underestimate of the LEZ effect. We 
cannot conclude, therefore, on the direction of the potential bias. 

The data analyzed in this study are the only available 
longitudinal measuring data to investigate the development of 
nitrogen oxide concentrations before and after the introduction of 
LEZs in Germany. We conclude that the material used can be 
considered as "data best available". Interpretations are limited, 
however, because spatial representativeness of the measuring sites 
can be disputed. It is unknown whether these data can be used to 
reliably estimate the exposures of citizens living in the LEZs. Since 
this is not only a problem of German measuring networks but an 
issue on the European level a research project was started to 
investigate the representativeness of measurement sites [78]. The 
authors concluded that measurements at the background stations 
are of greater importance than the data collected at the hot spots 
(traffic stations). Other limitations of hot spot data result from the 
fact that the citizens living in the LEZ area spend most of their 
time indoors and that indoor pollution data differ from hot spot 
outdoor concentrations [90,91]. 

Conclusions 

This is the first comprehensive approach to assess effects of 
LEZs on NO2, NO and NO^ concentrations with the help of 
measurement data on the Federal level in Germany. Reductions 
due to introducing LEZs of pollutant group 1 were estimated to be 
limited by 2 Hg/m"* (or 4%). The 4-week averages of NO2 
concentrations at index stations after LEZ introduction were found 
to be 55 |ig/m'' (median and mean) or 82 IXg/m"* (95th percentile). 
The NO2 concentration limit [24] enforced in Germany since 
2010 is 40 \xg/m'^ (1 year average). Concerning the expenditure of 
regulations and controls which are required to introduce and 
operate LEZs in cities, the proven impact of LEZs on the 
reduction of NOj ambient air concentrations with at a maximum 
of 4% in the first phase is very small. 

Supporting information 

File SI Contains Tables S1-S4 and Figures S1-S19. 
Table SI: Detailed results on NO2 - quadruplet analyses by linear 
(additive) log-linear (multiplicative) regression models. Table S2: 
NO: Log-linear (multiplicative) model 1 evaluating the quadru- 
plets of pooled continuous and diffiise sampler NO-measurements. 
Regression coefficient, robust standard errors of coefficient, t- 
statistic, t^\-o-sidcd P-\'alue, and 9,5 ""(i -confidence interval of 
coefficient. The relative LEZ effect estimate is given by the 
coefficient E (<1: concentration is lowered by LEZ). Table S3: 
NOx: Quadruplets of pooled continuous and diffuse sampler NOx- 
measurements: index stations (Ind), reference stations (Ref) before 
(pre) and after (post) introduction of LEZ. Ind.diff and Ref diff 
denote differences between index measurements and between 
reference measurements (negative post-pre differences indicate 



lower values after introduction of LEZ). Table S4: NO^: Linear 
(additive) model 1 evaluating the quadruplets of pooled continuous 
and diffuse sampler NO^-measurements. Regression coefficient, 
robust standard errors of coefficient, t-statistic, two-sided P-value, 
and 95%-confidence interval of coefficient. The absolute LEZ 
effect estimate is given by the coefficient E in |J.g/m'^ (<0: 
concentration is lowered by LEZ). Table S5: NO^: Log-linear 
(multiplicative) model 1 evaluating the quadruplets of pooled 
continuous and diffuse sampler NOx-measurcmcnts. Regression 
coefficient, robust standard errors of coefficient, t-statistic, two- 
sided P-value, and 9.5'X)-confidence interval of coefficient. The 
relative LEZ effect estimate is given by the coefficient E (<1: 
concentration is lowered by LEZ). Figure SI: Low emission zone 
Herrenberg (marked area), implemented in 2009-01-01 (mod- 
ified from www.map24.de). One index station: 1)DEBW135 
Hindenburger StraRe, no NO, no NOx. One reference station 
outside the low emission zone: 2)DEBW112 Gartringen (not 
included in the figure since located approx. 5 km north of low 
emission zone). Figure S2: Low emission zone Ilsfeld (marked 
area), implemented in 2008-03-01 (modified from www.map24. 
de). One index station: 1)DEBW133 Konig-Wilhelm-StraBe, no 
NO, no NOx. One reference station outside the low emission zone: 
2)DEBW034 Waibhngen (not included in the figure since located 
approx. 24 km south of low emission zone). Figure S3: Low 
emission zone Karlsruhe (marked area), implemented in 2009- 
01-01 (modified from www.map24.de). One index station: 
1)DEBW126 Kriegsstral?e, no NO2, no NOx. Two reference 
stations outside the low emission zone: 2)DEBW001 Kaiisruhe- 
Mitte 3)DEBW004 Eggenstein (not included in the figure since 
located approx. 6 km north of low emission zone). Figure S4: 
Low emission zone Ludwigsburg (marked area), implemented in 
2008-03-01 (modified from www.map24.de). Two index stations: 

1) DEBW024 Weimar-/Schweizerstr£iBe 2)DEBW017 FriedrichstraBe. 
One reference station outside the low emission zone: 3)DEBW034 
WaibHngen (not included in the figure since located approx. 7 km 
south east of low emission zone). Figure S5: Low emission zone 
Mannheim (marked area), implemented in 2008-03-01 (modified 
from www.map24.de). Two index stations: 1)DEBW006 Mannheim- 
Mitte 2)DEBW098 Friedrichsring U2. Two reference stations outside 
the low emission zone: 3)DEBW005 Mannlieiiu Nord (not included in 
the figure since located approx. 4 km north of low emission zone) 

4) DEBW007 Mannheim-Siid (not included in the figure since located 
approx. 5 km south of low emission zone). Figure S6: Low emission 
zone Reutlingen (marked area), implemented in 2008-03-01 
(modified from www.map24.de). One index station: 1)DEBW027 
EbertstraBe. Two reference stations outside the low emission zone: 

2) DEBW042 Bemhausen (not included in the figure since located 
approx. 20 km north of low emission zone) 3)DEBW117 Gartringen 
(not included in the %ure since located approx. 28 km west of low 
emission zoni-). Figure S7: Low emission zone Stuttgart (marked 
area), implemented in 2008-03-01 (modified from v™w.map24.de). Six 
index stations: 1)DEBW011 Zuffenhausen 2)DEBW()13 SeubcrstniBe 

3) DEBW099 Amulf-Klett-Platz 4)DEBW116 Hohenheimer StraBe 

5) DEBW118 Am Neckartor 6)DEBW134 Waiblinger StraBe. Two 
reference stations outside the low emission zone: 7)DEBW034 
Waiblingen 8)DEBW042 Bemhausen (not included in the figure since 
located approx. 2 km south of low emission zone). Figure S8: Low 
emission zone Tubingen (marked area), implemented in 2008-03-01 
(modified from www.map24.de). One index station: 1)DEBW107 
DerendingerstraBe. One reference station outside the low emission 
zone: 2)DEBW1 12 Gartringen (not included in the figure since located 
approx. 1.5 km north west of low emission zone). Figure S9: Low 
emission zone Augsburg (marked area), implemented in 2009-07-01 
(modified from www.map24.de). Three index stations: 1)DEBY007 
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Bourges-Hatz 2)DEBY110 KarlstraBe 3)DEBY006 Konigsplatz. One 
reference station outside the low emission zone: 4)DEBY099 LfU (not 
included in the figure since located approx. 3 km south of low emission 
zone). Figure SIO: Low emission zone Munich (marked area), 
implemented in 2008-10-01 (modified from www.map24.de). Five 
index stations: 1)DEBY037 Stachus 2)DEBY039 LothstraBe 
3)DEBY085 Luise-Kiesselbach-Platz 4)DEBY114 PrinzregentenstraBe 
e 5)DEBY1 15 Landshuter Allee. Three reference stations outside the 
low emission zone: 6)DEBY043 Moosach, no PM,,, 7)DEBY089Jo- 
hanneskirchen 8)DEBY1()9 Andechs/Rothenfeld (not included in the 
figure since located approx. 27 km south west of low emission zone). 
Figure SI la: Low emission zone Berlin Blume-Messnetz 
(marked area), implemented in 2008-01-01 (modified from www. 
map24.de).Five index stations: 1)DEBE018 B Schoneberg-Belziger 
StraBe 2)DEBE034 B Neukolln-NansenstraBe 3)DEBE064 B Neukolln- 
Karl-Marx-StraBe 76 4)DEBE065 B Friedrichshain-Frankfurter Allee 
5)DEBE067 B Hardenbergplatz. Nine reference stations outside the 
low emission zone: 6)DEBE061 B Steglitz-SchUdhomstraBe 

7) DEBE062 B Frohnau, Funkturm (not included in the figure since 
located approx. 1 3 km north of low emission zone) 8)DEBE063 B 
NeukoUn-SilbersteinstraBe) 9)DEBE066 B Karlshorst-RheingoldstraBe, 
no PMio (not included in the figure since located approx. 5 km east of 
low emisskm zone) 10)DEBE010 B Wedding-Amrumer StraBe 

1 1) DEBE027 B Marienfelde-Schichauweg (not included in the figure 
since located approx. 8 km south of low emission zone) 1 2)DEBE032 B 
Grunewald (not included in the figure since located approx. 4 km 
south west of low emission zone) 1 3)DEBE051 B Buch (not included in 
the figure since located approx. 1 2 km north east of low emission zone) 

14) DEBE()56 B Friedrichshagen (not included in the figure since 
located approx. 14 km south east of low emission zone). Figure SI lb: 
Low emission zone Berlin RUBIS-Messnetz (marked area), 
implemented in 2008-01-01 (modified from www.map24.de). Ten 
index stations: 1)DEBE530 HauptstraBe 30 2)DEBE504 BeusselstraBe 
66 3)DEBE537 Alt Moabit 63 4)DEBE545 Sonnenallee 68 
5)DEBE547 Landsberger AUee 6-8 6)DEBE517 Neukolln- 
NansenstraBe 7)DEBE519 Friedrichshain-Frankfiarter Allee 

8) DEBE555 Herrmannplatz Lateme 21 9)DEBE562 FriedrichstraBe 
Lateme 156 10)DEBE525 Leipziger StraBe 32. Twelve reference 
stations outside the low emission zone: 1 1) DEBE501 Berliner Allee 118 

12) DEBE577 Buch, no NO, no NO^ (not included in the figure since 
located approx. 12 km north of low emission zone) 13)DEBE507 
Griinauer StraBe 4 (not included in the figure since located approx. 
9 km south east of low emission zone) 14)DEBE539 SchloBstraBe 29 

15) DEBE542 Tempelhofer Damm 148 16)DEBE513 SpreestraBe 2 
(not included in the figure since located approx. 5 km south east of low 
emission zone) 1 7)DEBE5 14 Alt Friedrichsfelde 8a (not included in the 
figure since located approx. 3 km east of low emission zone) 
18)DEBE521 StegHtz-SchildhomstraBe 19)DEBE559 Buschkrugallee 
Lateme 3 20)DEBE522 Neukolln-SilbersteinstraBel 21)DEBE573 
BadstraBe 22)DEBE576 Spandau, KlosterstraBe 12 (not included in 
the figure since kx'at(;d approx. 6 km west of low emission zone). 
Figure S12: Low emission zone Frankfurt a.M. (marked area), 
implemented in 2008-10-01 (modified from www.map24.de). One 
index station: 1)DEHE041 Frankfurt-Friedb.Ldstr. Three reference 
stations outside the low emission zone: 2)DEHE008 Frankfurt-Ost 
3)DEHE01 1 Hanau (not included in the figure since located approx. 
13 km east of low emission zone) 4)DEHE005 Frankfurt-Hochst. 
Figure S13: Low emission zone Hannover (marked area), 
implemented in 2008-01-01 (modified from www.map24.de). One 
index station: 1)DENI048 Hannover Verkehr. Four reference stations 
outside the low emission zone: 2)DENI054 Hannover 3)DENI011 
Braunschweig, Broizemer Steinberg (not included in the figure since 
located approx. 49 km east of low emission zone) 4)DENI041 
Weserbei^land/Rinteki, Brugfeldsweide (not included in the figure 



since located approx. 48 km south west of low emission zone) 

5) DENI052 AllertalAValsrode, Auf dem Kamp 8 (not included in tiie 
figure since located approx. 47 km north of low emission zone). 
Figure S14: Low emission zone Dortmund (marked area), 
implemented in 2008-10-01, but Brackelerstr. 2008-01-01 (modified 
from www.map24.de). Four index stations: 1)DEN\V1()1 SteinstrtiBe 

2) DENW136 Brackeler StraBe 3)DENW184 Westialendamm 190, no 
NO, no NOx, no PMio 4)DENW185 Rheinlanddamm 5-7, no NO, 
no NOx, no PMk). Four reference stations outside the low emission 
zone: 5)DEN\V0()2 Datteln-Hagem (not included in the figure since 
located approx. 15 km iiortli ^\ (:st of low emission zone) 6)DEN\V008 
Do-Eving 7)DENW029 Hattingen, An der Becke (not included in the 
figure since located approx. 19 km soutii west of low emission zone) 

8) DIiNW179 Schwerte (not included in the figure since located 
approx. 8 km south of low emission zone). Figure S15: Low emission 
zone Duisburg (marked area), implemented in 2008-10-01 (modified 
from www.map24.de).Three index stations: 1)DENV\''034 Duisburg- 
Walsum 2)DENW040 Dukburg-Buchholz 3)DENW112 Kardinal- 
Galen-StraBe. One reference station outside the low emission zone: 
4)DENW038 45476 Mtihlheim, NeustadtstraBe (not included in the 
figure since located approx. 5 km east of low emission zone). Figure 
S16: Low emission zone Diisseldorf (marked area), implemented in 
2009-02-15 (modified from www.map24.de). Two index stations: 

I) DENW082 ComeHusstraBe 2)DENW216 Diisseldorf-Bilk, no NO, 
no NOx, no PMig. Four reference stations outside the low emission 
zone: 3)DENW042 Krefeld-Linn (not included in the figure since 
located approx. 14 km north west of low emission zone) 4)DENW071 
Dilsseldorf-Lcirkk (not inc:luded in die figure since located approx. 
3 km west of low emission zone) 5)DEN'\¥078 Ratingen-Tiefenbroich 
(not included in the figure since kicated approx. 6 km north east of low 
emission zone) 6)DENW1 16 Krefeld Hafen (not included in the figure 
since located approx. 12 km north west of low emission zone). Figure 
S17: Low emission zone Essen (marked area), implemented in 2008- 
10-01 (modified from www.map24.de). Eight index stations: 1)DEN- 
W043Ost Steeler StraBe 2)DENW134 Gladbecker StniBe 

3) DENW135 Hombrucher StraBe 4)DENW161 AlfredstraBe 9/11, 
no NO, no NO,, no PMio 5)DENW168 Gladbecker StraBe 245, no 
NO, no NO,, no PMm 6)DENW169 In der Baumschule 7, no NO, no 
NO,, no PMio 7)DENW171 HombrucherstraBe 21/23, no NO, no 
NO,, no PMio 8)DEN\V215 HausackerstraBe 1 1, no NO, no NO,, no 
PMio. Three reference stations outside the low emission zone: 

9) DENW024 Essen-Vogelheim 10)DENW029 Hattingen-Blanken- 
stein (not included in the figure since located approx. 10 km south east 
of low emission zone), 1 1) DENW162 BrilckstraBe 29, no NO, no NO,„ 
no PMio (not included in the figure since located approx. 4 km south of 
low emission zone). Figure S18: Low emission zone Cologne 
(marked area), implemented in 2008-01-01 (modified from www. 
map24.de). Seven index stations: 1)DENW148 JustinianstraBe 13-15, 
no NO, no NO,, no PMjo 2)DENW151 Neumarkt 25, no NO, no 
NO,, no PMio 3)DEN\V153 TunisstraBe/Elstergasse, no NO, no 
NO,, no PMio 4)DENVV164 Hohenstaufenring 57A, no NO, no NO,, 
no PMio 5)DEN\V198 Gereonsdiiesch 21, no NO, no NO,, no PMk, 

6) DENW211 Clevischer Ring 3 7)DENW212 Turiner StraBe. Four 
reference stations outside the low emission zone: 8)DENW053 
Cologne-Chorweiler (not included in the figure since located approx. 
9 km north west of low emission zone) 9)DENW058 Hiirth (not 
included in the figure sin[:e located approx. 7 km south west of low 
emission zone), 1 ())DEN\V059 Cologne-Rodenkirchen (not included in 
the figure sinc:e kK:ated approx. 4 km south of low emission zone), 

I I) DENW079 Leverkusen-Manfort (not included in the figure since 
located approx. 7 km north of low emission zone). Figure S19: Low 
emission zone Wuppertal (marked area), implemented in 2009-02- 1 5 
(modified from www.map24.de). Two index stations: 1)DENW114 
Wuppertal-Langerfeld, no NO, no NO, 2)DENW189 Wuppertal 
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Gatlie. Two reference stations outside the low emission zone: 
3)DENW029 Hattingen-Blankenstein (not included in the figure since 
located approx. 13 km north of low emission zone) 4)DENW080 
Solingen-Wald (not included in the figure since located approx. 5 km 
south west of low emission zone). 
(DOCX) 

Acknowledgments 

We thank the project steering committee and the following stale authorities 
for the intense and helpful discussion and the supply of the raw data. In 
specilic: Peter Bruckniann and Reinhold Beier, Landesamt fiir Natur-, 
Umwelt- und Verbraucherschutz NRW, Abteilung 4, Luftqualitat, 
Gerauselie, Erschiitterungen, Strahlensehutz; Stefan Jacobi and Wieslawa 
Stec-Lazaj, Hessisches Landesamt fiir Umwelt und Geologic (HLUG), 
Abteilung I Immissions- und Strahlensehutz; Michael Koster and Andreas 
Hainsch, Staadiches Gewerbeaufsichtsamt HUdesheim, Abteilung 4 - 

References 

1. LEZEN (2013) Low emission zone in Europe network. Available from: http:// 
w^vAv.lowcmissionzoncs.cu. Accessed 2014July 8. 

2. EC (2007) Regulation (EC.'.) No 71.V2007on type approval of motor vehicles 
with respect to emissions from light passenger and commercial vehicles (Euro 5 
and Euro 6) and on access to vehicle repair and maintenance information. 1—16. 
Available from: http: / / cur-lex. europ a. eu/ Lex UriServ/ LcxUriServ. 
do?uri-OJ:L:2007:171:0001:0001:EN:PDF. Accessed 2014July 8. 

3. Umweltbundesamt (2014) Umweltzonen in Deutschiand. Available from: 
http:/ / www.umweltbundesamt.de/themen/luft/luftschadstofFe/ feinstaub/ 
umweltzonen-in-deutsehland. Accessed 20 14 July 8. 

4. Colvile RN, Hutchinson EJ, Mindell JS. Warren RF (2001) The transport sector 
as a source of air pollution. Atmospheric Environment 35: 1537-1565. 

5. Bruckmann P, Eutz M (2010) V'crbessern Umweltzonen die Luftqualitat? In: 
Verband der Automobilindustrie (VDA), editor. 12. Technischer Kongress, 24. 
und 25. Miirz. Forum am Schlosspark, Ludwdgsburg: Henrich Druck+Medien 
GmbH. pp. 299-311. 

6. Bruckmann P, Wurzler S, Brandt A, Vogt K (2011) Erfahrungen mit 
Umweltzonen in Nordrhein-Westfalen. In: Bundesamt fur Strahlensehutz, 
Bundesinstitut fur Risikobewertung, Robert Koch-Institut, Umweltbundesamt, 
editor. UMID Umwelt und Mensch - Informationsdienst. Berlin, pp. 27—33. 

7. Kaesoh L (2011) Umweltzonen in Europa und in Deutschiand. In: Bundesamt 
fiir Sirahlrnschiilz, Biindrsinsliliii fiir Risikohewcrtiing, Robert Koch-Institut, 
Umweltbundesamt, editor. UMID Umwelt und Mensch - Informationsdienst. 
Berlin, pp. 5-10. 

8. Morfeld P, Keil U, Spallck M (2013) The European "Year of the Air": fact, fake 
or vision? Archives of Toxicology 87: 2051—2055. 

9. Turner MC, Krewski D, Pope CA, Chen Y, Gapstur SM, et al. (2011) Long- 
term ambient fine particulate matter air pollution and lung cEincer in a large 
cohort of never-smokcrs. American Journal of Respiratory and Critical Care 
Mcdicme 184: 1374-1381. 

10. Hock G, Krishnan RM, Beclen R, Peters A, Ostro B, et al. (2013) Long-term air 
pollution exposure and cardio- respiratory mortality: a review. Environmental 
Health 12: 1-15. 

1 1 . Hystad P, Demers PA, Johnson KC, Garpiano RM, Brauer M (20 1 3) Long-term 
residential exposure to air pollution and lung cancer risk. Epidemiology 24: 762- 
772. 

12. Raaschou-Nielsen O, Andersen ZJ, Beelen R, Samoli E, Stafoggia M, et sd. 
(2013) Air pollution and lung cancer incidence in 17 European cohorts: 
prospective analyses from the European Study of Cohorts for Air Pollution 
Effects (ESCAPE). The Lancet Oncology 14: 813-822. 

13. Schuster G (2009) Umweltzonen: Nutzen weitcrhin umstritten. Dtsch Arztebl 
106: A228. 

14. Eikmann T, Herr G (2009) 1st die Einfuhrung von Umweltzonen tatsachlich eine 
sinnvolle Mafinahme zum Schutz der Gesundheit der Bevolkerung? Umweltmed 
Forsch Prax 14: 125-126. 

15. Friedrich B (2008) Umweltzonen. Strafienverkehrstechnik 11: 673. 

16. Zellner R, Kuhlbusch TAJ, Diegmann V, Herrmann H, Kasper M, et al. (2009) 
Feinstaubc und Umweltzonen. Available from: Available from: www.processnet. 

org/ dechema_media/ Downloads /Positionspapiere/Stellungnahme 

Feinstaeube.pdf. Accessed 2014 July 8. 

17. Downs SH, Schindler C, Liu LJ,'Kcidel D, Bayer- Ogle. shy L, et al. (2007) 
Reduced exposure to PMIO and attenuated age-related decline in lung function. 
The New England Journal of Medicine 357: 2338-2347. 

18. Gauderman WJ, Avol E, Gilliland F, Vora H, Thomas D, et al. (2004) The effect 
of air pollution on lung development from 10 to 18 years of age. The New 
England Journal of Medicine 351: 1057-1067. 

19. Lim SS, Vos T, Flaxman AD, Danaei G, Shibuya K, et al. (2013) A comparative 
risk assessment of burden of disease and injury attributable to 67 risk factors and 
risk factor clusters in 21 regions, 1990-2010: A systematic analysis for the Global 
Burden of Disease Study 2010. The Lancet 380: 2224-2260. 



Zentrale UnlcrstuLzungsstelle fur Luftrcinhallung und (jicfahrstofrc (ZUS 
LG); Martin Lutz and Arnold KciLschau, Berliner Senatsverwaltung fur 
Gesundheit, Umwelt und Verbrauchcrsehulz, Abteilung III Umweltpolitik, 
Referat III D Immissions schutz; Heinz Ott, Bayerisches Landesamt fur 
Umwelt Abteilung 2, Referat 2, 4 Luftgiitemessungcn Siidbayern, 
Luftreinhaltung beim Verkehr; Werner Scholz and Ghristiane Lutz- 
Holzhauer, Landesanstalt fiir Umwelt, Messungen und Naturschutz 
Baden-Wiirttemberg (LUBW), Referat 33 Luftqualitat; Ralf Wehrse and 
Jan Osmers, Bremer Senator fiir Umwelt, Bau, Verkehr und Europa 
Fachbereich Umwelt, Abteilung 2, Umweltwirtschafts-, Klima- und 
Ressourcenschutz, Referat 22 Immissionsschutz. 

Author Contributions 

Analyzed the data: PM. Contributed reagents/materials/analysis tools: 
PM. Contributed to the writing of the manuscript: PM DO MS. Data 
Organization; MS DG PM. Health Aspects: MS DG. 



20. Nawrot TS, Perez L, Kunzli N, Munters E, Nemery B (2011) Public health 
importance of triggers of myocardial infarction: a comparative risk assessment. 
The Lancet 377: 732-740. 

21. Samct JM, Dominici F, Curriero FC, Goursac 1, Zeger SL (2000) Fine 
particulate air pollution and mortality in 20 U.S. cities, 1987-1994. The New 
England Journal of Medicine 343: 1742-1749. 

22. Latza U, Gerdes S, Baur X (2009) Effects of nitrogen dioxide on human health: 
systematic review of experimental and epidemiological studies conducted 
between 2002 and 2006. International Journal of Hygiene and Environmental 
Health 212: 271-287. 

23. EC (2008) Directive 2008/50/EC of the European Parliament and of the 
councUof 21 May 2008 on ambient air quality and cleaner air for Europe. 1—44. 
Available from: http:/ / cur-lex. europa. eu/LcxUriServ/LexUriServ. 
do?uri^GJ:L:2008:152:0001:0044:E\:PDF. Accessed 2014 July 8. 

24. EC (1999) (Joiineil Directi\'e 1999/30/ECJ of 22 April 1999 relating- lo limit 
values for sulphur dioxide, nitrogen dioxide and oxides of nitrogen, particulate 
matter and lead in ambient air. 41-60. Available from: http://eur-lex. europa. 
eu/LexUriServ/LexUriSer\-.do?uri ^ OJ:L:1999:163:0041:0060:EN:PDF. Ac- 
cessed 2014 July 8. 

25. Umweltbundesamt (2012) Luftqualitat 2011 - Feinstaubepisoden pragten das 
Bild. Dessau- RoGlau. Available from: www.umweltbundesamt.de/uba-info- 
medien/421Lhtml. Accessed 2014 July 8. 

26. Giannouli M, Kalognomou E-A, Mellios G, Moussiopoulos N, Samaras Z, et al. 
(2011) Impact of European emission control strategies on urban and local air 
quality. Atmospheric Environment 45: 4753-4762. 

27. European Environment Agency (2014) Indicator: Excredance 492 of air equality 
limit values in urban areas. Copenhagen, Denmark: European Environment 
Agency. Available from: http:/ /wwvv.eea.europa.eu/ data-and-maps/indicators/ 
exceedancc-of-air-quality-limit- 1 / exceedancc-of-air-quality-limit-5. Accessed 
20 14 July 16. 

28. European Environment Agency (2014) AirBase - The European air Quality 
database. Copenhagen, Denmark. Available fr-om: http://acm.eionet.europa. 
eu/databases/airbase. Accessed 2014 July 8 http://www.eea.europa.eu/data- 
and-maps/ data/airbase-the-european-air-quality-database-8. Accessed 20 1 4 Ju- 
ly 8. 

29. Rauterberg-Wulff A, Lutz M (2011) Wirkungsuntersuchungen zur LImwcltzone 
Berlin. In: Bundesamt fiir Strahlensehutz, Bundesinstitut fiir Risikobewertung, 
Robert Koch-Institut, Umweltbundesamt, editor. UMID Umwelt und Mensch - 
Informationsdienst. Berlin, pp. 11—18. 

30. Lutz M, Rauterberg-Wulff A (2010) Berlin's low emission zone - top or flop? 
14th ETH Conference on Combusion (Generated Particles. 

31. Lutz M, Rauterberg-Wulff A (2009) Ein Jahr Umweltzone Berlin: Wirkungsun- 
tersuchungen. Berlin. 1—3 1 . Available from: www.berlui.de/ sen/ umwelt/ 
luftqualitaet/de/luftreuihalteplan/download/umweltzone_ ljahr_bericht.pdf. 
Accessed 20 14 July 8. 

32. Morfeld P, Groneberg D, Spallck M (2014) Effectiveness of low emission zones: 
Analysis of the changes in fine dust concentrations (PM]o) in 19 (ierman cities. 
Pneumologie 68: 173-186. 

33. Vogt R, Kessler C, Schneider C (2010) Stadtischc NOj Luftqualitat: 
Quelle nanalyse und zukiinftige Entwicklung. In: Verband der Automobilindus- 
trie (VDA), editor. 12 Technischer Kongress 2010, 24 und 25 Marz. Forum am 
Schlosspark, Ludwigsburg: Henrich Druck-I-Medien GmbH. pp. 287-297. 

34. US EPA (2008) Integrated science assessment for oxides of nitrogen - health 
criteria (fmal report). EPA/600/R-08/071. 260. Available fi-om: http://www. 
epa.gov/ord/htm/whatsnew.htm. Accessed 2014July 8. 

35. EC (2001) Directive 2001/81/EC of the European Parliament and of the 
council of 23 October 2001 on national emission ceilings for certain atmospheric 
pollutants. 22-30. Available from: http:/ / eur-lex. europa. eu/LexUriServ/ 
LexUriServ.do?uri = OJ:L:2001:309:0022:0030:EN:PDF. Accessed 2014 July 8. 



PLOS ONE I www.plosone.org 



16 



August 2014 I Volume 9 | Issue 8 | e102999 



Effectiveness of LEZs 



36. Atkinson R\V, Barratt B, Armstrong B, Anderson HR, Bccvcrs SD, ct al. (2009) 
The impact of the congestion charging scheme on ambient air pollution 
concentrations in London. Atmospheric Environment 43: 5493—5500. 

37. Bruckmann P, Brandt A, Wurzler S, Vogt K (201 1) Verbessem Umweltzonen 
die Luftqualitat? Neue Entwiddungen bei der Messung und Beurteilung der 
Luftqualitaet Fachtagung der Kommission Reinhaltung der Lufl im VDI und 
DIN-Normcnausschuss KRdL, Badcn-Badcn, 11-12 Mai. pp. 3-24. 

38. Umwcltbundrsamt (2014) Stationsdatcnbank dcs Umwrltbundrsamtes. Avail- 
able from: http;/ /www.cnv-it.dc/ stationcn/pubtic/ stationList.do. Accessed 
20 14 July 8. 

39. Pfclicr U (2010) Mcsstcchnik fur Stickstolfdioxid (NO^). In: Kommission 
Reinhaltung der Luft im VDI und DIN-Normcnausschuss KRdL, editor. 
Stickstoffdioxid und Partikel (PMa^s/PMio). pp. 1 13-122. 

40. Umweltbundesamt (2014) Luft und LuflreinhEiltung. Available from: http:// 
www.umweltbundesamt.de/ themen/luft/ messenbeobachtenueberwachen/ 
messgeraete-messverfahren. Accessed 2014 July 8. 

41. Pfeffer U, Bier R, Zang T (2006) Measurements of nitrogen dioxide with 
diffusive samplers at traffic-related sites in North Rhine -Westphalia (Germany). 
Gcfahrstoffc, Reinhaltung der Luft 1/2: 38-44. 

42. Pfeffer U, Zang T, Rumpf E-M, Zang S (2010) Calibration of diffusive samplers 
for nitrogen dioxide using the reference method - Evaluation of measurement 
uncertainty. Gefahrstoffe, Reinhaltung der Lui't 11/12: 500-506. 

43. DIN EN 13528-3 Aufienluftqualitat (2004) Passivsammler zur Bestimmung der 
Konzentrationen von Gasen und Dampfen - Tail 3: Anleitung zur Auswahlt, 
Anwendung und Handhabung. Berlin: Beuth. 

44. Rothman KJ, Greenland S, Lash TL (2008) Modem Epidemiology. 3. ed. 
Philadelphia: Lippincott Williams & Wilkins. 

45. Morfeld P, Spallek M, Gronebcrg D (201 1) Zur Wirksamkcit von Umweltzonen: 
Design einer Studie zur Ermittlung der Schadstolfkonzcntrationsanderung fiir 
Staubpartikel (PMio) und andere ()ro|3('n dureh Einfiihrung von Umweltzonen 
in 20 deutschen Stadten. Zentralblatt fur Arbeitsmcdizin, Arbcitsschutz und 
Ergonomie 61: 148-165. 

46. Kelly F, Armstrong B, Atkinson R, Anderson HR, Barratt B, et al. (2011) The 
London low emission zone baseline study. Health Effects Institut 163. 1-96. 
Available from: http://pubs.healtheffects.org/vicw.phppid = 366. Accessed 2014 
July 8 

47. Allison PD (2009) Fixed effects regression models. Los Angeles: SAGE. 

48. Royston P, Sauerbrei W (2008) Multivariable model-buUdung. Chichester, 
England: John Wiley & Sons Inc. 1—303 p. 

49. Lenschow P, Abraham HJ, Kutzner K, Lutz M, PreuRJD, et al. (2001) Some 
ideas about the sorces of PMIO. Atmospheric Environment 35: S23-S33. 

50. Hock G, Meliefste K, Cyrys J, Lewne M, Bellander T, et al. (2002) Spatial 
variability of fine particle concentrations in three European areas. Atmospheric 
Environment 36: 4077-4088. 

51. Senn S (1997) Editorial: regression to the mean. Statisticsd Methods in Medical 
Research 6: 99-102. 

52. Altman DG, Bland AE (1983) Measurement in medicine: the analysis of method 
comparison studies. The Statistician 32: 307-317. 

53. Bland JM, Altman DG (1986) Statistical methods for assessing agreement 
between two methods of clinical measurement. The Lancet 1: 307—310. 

54. GillJS, Zezulka AV, Bccvcrs DG, Davies P (1985) Relation between initial blood 
pressure and its fall with treatment. The Lancet 1: 567 569. 

55. Bland JM, Altman DG (1994) Regression towards the mean. British Medical 
Journal 308: 1499. 

56. Bland JM, Altman DG (1994) Some examples of regression towards the mean. 
British Medical Journal 309: 780. 

57. Stigler SM (1997) Regression towards the mean, historically considered. 
Statistical Methods in Medical Research 6: 103-114. 

58. Bamett AG, van der Pols JC, Dobson AJ (2005) Regression to the mean: what it 
is and how to deal with it. International Journal of Epidemiology 34: 215—220. 

59. Twisk JVVR (2004) Applied longitudinal data analysis for epidemiology. 
Cambridge: Cambridge University Press. 62-77 p. 

60. Klingncr M, Sahn E, Anke K, Hoist T, Rost J, ct al. (2006) Rcduktionspo- 
tenziale vcrkchrsbcschrankcndcr MaBnahmcn in Bezug zu meteorologisch 
bedingten Schwankungen der PMiq- und NOx-Immissionen. Gefahrstoffe, 
Reinhaltung der Luft 66: 326-334. 

6 1 . Builtjes P, Jdifl W, Theloke J, Thiruchittampalam B, van der Gon HD, et al. 
(2012) Strategien zur Verminderung der Feinstaubbelastung. 1—160. Available 
from: http:/ / www.umweltbundcsamt.de/ sites /default /files/ medien/ 461/ 
publikationen/ 1268.pdf Accessed 2014 July 8 

62. Reimer E, Seherer B (1992) An Operational Meteorological Diagnostic Svstem 
for Regional Air Pollution Analysis and Long Term Modeling. In: Dop H, 
Kallos G, editors. Air Pollution Modeling and Its Application IX: Springer US. 
pp. 565-572. 

63. Stern R (2013) Anwendung des REM-CALGRID-Modells auf die Ballungs- 
raiune Berlin, Miinchen und Ruhrgebiet. Berlin: Freie Universitat Berlin, 
Institut flir Meteorologie. Tropospharische Umweltforschung. 67/2013. 1-95. 
Available from: http://www.umweltbundesamt.de/sites/default/files/medien/ 
46 1 / publikationen/texte_67_20 1 3_appelhans_m 1 4_komplett_0.pdf. Accessed 
2014 July 8. 

64. Kersehbaumer A, Reimer E (2003) Erstcllung der Mcteorologischcn Eingangs- 
daten fiir das REM/Calgrid-Modell: ModeUregion Berlin-Brandenburg. Ab- 
schlussbericht zum UBA-Forschungsvorhaben 29943246. Freie Universitat 



Berlin, Institut fur Meteorologie. Available from: http://opus.kobv.de/zlb/ 
volltexte/2009/7447/pdiy3729.pdi: Accessed 2014July 17. 

65. Kersehbaumer A (2010) Abhangigkeit der RCG-Simulationen von unterschie- 
cUichen meteorologischen Treibern. Forschungs-Teilbericht an das Umweltbun- 
desamt, im Rahmen des PAREST-Vorhabens: FKZ 206 43 200/1 "Strategien 
zur Verminderung der Feinstaubbelastung". Institut fur Meteorologie der Freien 
Universitat Berlin. Available from: http: / / www.iunweltbundesEunt.de/ sites/ 
default /filcs/mcdien/ 46 1 / publikationen/ tcxte_55_201 3_appelhans_m02_ 
komplctt_0_0.pdf Accessed 2014 July 17. 

66. Graedcl TE, Crutzen PJ (1994) Chemie der Atmosphare. Heidelberg: Spektrum 
Akademischer Verlag. 

67. Morfeld P, Stem R, Builtjes P, Groneberg DA, Spallek M (2013) Einrichtung 
einer Umweltzone und ihre Wirksamkeit auf die PMiQ-Feinstaubkonzentration - 
eine PUotanalyse am Beispiel Miinchen. Zentralblatt fiir Arbeitsmcdizin, 
Arbcitsschutz und Ergonomie 63: 104—115. 

68. Stern R (2003) Entwicklung und Anwendung des chemischen TransportmodcUs 
REM/CALGRID. Abschlussbcricht zum Forschungs- und Entwicklungsvorha- 
ben 298 41 252 des Umweltbundesamts "Modellierung und Prtifung von 
Strategien zur Verminderung der Belastung dureh Ozon". Available from: 
http:/ /www. umweltbundesamt.de/ sites/ default/files/medicn/publikation/ 
long/3604.pdf Accessed 2014July 17. 

69. Stem R (2004) GroBraumige PMio-AusbrcitungsmodeUiemng: Abschatzung der 
gcgenwartigen Immissionsbelastung in Europa imd Prognose bis 2010. KRdL- 
Experten-Fomm "Staub und Staubinhaltsstoffe", 2004-11-11/10. VDI-KRdL- 
Schriftenrcihe 33. 

70. Stern R (2004) Weitere Entwicklung und Anwendung des ehemisehcn 
TransportmodcUs REM-CALGRID fiir die bundeseinhcitliche Umsctzung der 
EU-Rahmenriehtlinie Luftqualitat und ihrcr Tochterriehtlinicn. Abschlussbcr- 
icht zum FuE-Vorhaben 201 43 250 des Umweltbundesamts "Anwendung 
modeUgestutzter BeurteUimgssysteme fiir die bundeseinheithche Umsetzung der 
EU-Rahmenrichdinie LuftquaHtat und ihrer Tochterrichtlinien". Available 
from: http:/ /www. umweltbundesamt.de/ sites/ default/ files/ medien/ 
publikation/long/36I0.pdf. Accessed 2014Julyl7. 

71. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical 
and powerful approach to multiple testing. Journal of the Royal Statistical 
Society, Scries B 57: 289-300. 

72. StataCorji (2009) Stata: Release II. Statistical Software. College Station, TX: 
StataCorp LP. 

73. Vickers AJ, Altman DG (2001) Statistics notes: Analysing controlled trials with 
baseline and follow up measurements. BMJ 323: 1123-1124. 

74. Neter J, Wasserman W, Kutner MH (1985) Applied linear statistical models. 
Regression, Analysis of variance and experimental designs. Homewood, Illinois: 
Richard Dr. Irwin. 

75. Puis T, Jager-Ambrozewiez (2012) Die Auswirkungen der Umweltzone in 
Frankfurt auf die N(.)2-lmmissionen. Koln: Institut der deutschen Wirtschaft 
Kdln. pp. 1-18. 

76. ZUS LG (2010) Bcwertung der Auswirkungen der Umweltzone Hannover auf 
Basis von Messdaten. HUdesheim: Zentrale Unterstiitzungsstelle Luftreinhaltung 
und Gefahrstoffe - Dezernat 42. 15. Available from: http://www.umwelt. 
niedersachsen.de/download/48880. Accessed 2014July 8. 

77. Umweltbundesamt (2014) Mediendatenbank. Available from: http://www. 
umweltbundesamt.dc/ publikationcn?kcys = PAREST&topic = AU&series = AIl&sort_ 
beflcombinc — field_date_monthly_value-l-DESC. Accessed 2014July 17. 

78. Duyzer J, van den Hout D, Zandveld P, van Ratingen S (2013) Representa- 
tiveness of air quality monitoring stations. TNO Research Report 2013 Rl 1055. 
TNO Utrecht (in print). 

79. Tonne C, Beevers S, Armstrong B, Kelly F, Wilkinson P (2008) Air pollution and 
mortality benefits of the London Congestion Charge: spatial and socioeconomic 
inequalities. Occupational and Environmental Medicine 65: 620-627. 

80. Briggs DJ (2008) A framework for integrated environmental health impact 
assessment of systemic risks. Environmental Health 7: 61-77. 

81. Cesaroni G, Boogaard H, Jonkcrs S, Porta D, Badaloni C, et al. (2012) Health 
benefits of trallie-related air pollution reduction in different socioeconomic 
groups: the effect of low-emission zoning in Rome. Occupational and 
Environmental Medicine 69: 133-139. 

82. Johansson C, Burman L, Forsberg B (2009) The effects of congestions tax on air 
quality and health. Atmospheric Environment 43: 4843^854. 

83. Boogaard H, Janssen NAH, Fischer PH, Kos GPA, Weijers EP, et al. (2012) 
Impact of low emission zones and local traffic policies on Eunbient air pollution 
concentrations. The Science of the Total Environment 435—436: 132—140. 

84. EU (2013) Low emission zones in Europe. Available from: http://www. 
lowemissionzones.cu/. Accessed 2014 July 8. 

85. HBEFA (2010) Handbook emission factors for road transport (HBEFA 3.1). 
Available from: http://www.hbefa.net/e/index.html. Accessed 2014July 8. 

86. Carslaw DC, Bccvcrs SD, Tate JE, Westmoreland EJ, Williams ML (2011) 
Recent evidence concerning higher NO^ emissions from passenger cars and light 
duty vehicles. Atmospheric Environment 45: 7053-7063. 

87. Beevers SD, Westmoreland E, de Jong MC, Williams ML, Carslaw DC (2012) 
Trends in NO^ and NO2 emissions from road traffic in Great Britain. 
Atmospheric Environment 54: 107-116. 

88. LANUV, Landcsamt fiir Natur, Umwclt und Vcrbraucherschutz NRW (2008) 
Vorgehensweise des LANUV zur Korrektur kontinuierlichcr PMIO-Mcssdaten 
im Luftmessnetz von NRW. Available from: http://www.lanuv.nrw.de/luft/ 
imniissionen/ber_trend/erlaeutemngen.pdf. Accessed 20 14 July 8. 



PLCS ONE I www.plosone.org 



17 



August 2014 I Volume 9 | Issue 8 | e102999 



Effectiveness of LEZs 



89. LUBW, Landcsanstalt fiir Umwclt, Mcssungcn und Naturschutz Badcn- 
Wurttemberg (2009) Untersuchung von masscnrelcvanten InhaltsstolFen in 
Feinstaub PMio- Karlsruhe. 1-72. Available from: www.lubw.baden- 
wuerttemberg.de/servlet/is/207409/untersuchung_massenrelevanten_ 
inh£iltsstoffen_feinstaub_pml0.pdf?command = downloadContent&filename = 
untersuchung_massenrelevanten_inhaltsstofren_feinstaub_pm 1 0.pdf. Accessed 
2014 July 8. 



90. Dons E, Int Panis L, Van Poppcl M, ThcunisJ, Willcms H, ct al. (201 1) Impact 
of time-activity patterns on personal exposure to black carbon. Atmospheric 
Environment 45: 3594-3602. 

91. Fischer PH, Hoek G, van Reeuwijk H, Briggs DJ, Lebret E, et al. (2000) Traffic- 
related differences in outdoor and indoor concentrations of particles and volatile 
organic compounds in Amsterdam. Atmospheric Environment 34: 3713-3722. 



PLOS ONE I www.plosone.org 



18 



August 2014 I Volume 9 | Issue 8 | e102999 



